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Preface 


The  purpose  of  this  study  was  to  determine  if  Box-Jenkins 
time  series  techniques  forecast  C-135  aircraft  recoverable 
spares  accurately.  Emphasis  has  been  placed  on  finding  the 
best  methods  for  forecasting  aircraft  recoverable  spares 
since  they've  been  greatly  underestimated  in  the  last  few 
years.  This  research  consisted  of  developing  Box-Jenkins 
transfer  function  models  which  used  past  flying  hours  to 
forecast  future  demand,  Box-Jenkins  univariate  models  which 
used  past  demand  and  its  relationship  over  time  to  forecast 
future  demand,  and  simple  exponential  smoothing  models  which 
emphasize  recent  demand  data  more  heavily  in  forecasting 
future  demand.  All  three  types  of  models  produced  very 
accurate  results  when  forecasting  demand  one  quarter  into 
the  future.  The  simple  exponential  smoothing  models  had 
the  most  accurate  results  of  the  three  for  the  quarter. 

This  thesis  would  not  have  been  possible  without  the 
assistance  of  certain  individuals.  A  special  thanks  is  owed 
to  Lieutenant  Colonel  Bruce  Christensen  whose  help  and 
guidance  as  my  advisor  was  greatly  appreciated.  I  would  also 
like  to  thank  Captain  Terry  Tong  and  Larry  Collins  for  helping 
me  obtain  flying  hours  and  demand  data.  Finally,  a  heartfelt 
thanks  goes  to  my  family  and  friends  for  their  love  and 
support  throughout  this  research  endeavor. 

Tammy  M.  Haight 
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Abstract 

In  recent  years.  Air  Staff  directed  a  comprehensive 
review  of  the  spares  forecasting  process  due  to  significant 
underestimates  of  spares  requirements.  The  purpose  of  this 
study  was  to  determine  if  time  series  forecasting  models 
could  accurately  forecast  demand  for  aircraft  recoverable 
spares.  Box-Jenkins  time  series  analysis  was  used  to  analyze 
and  develop  forecasting  models  for  ten  C-135  spares. 

Two  different  Box-Jenkins  models  were  developed  for  each 
spare  to  forecast  demand.  These  forecasts  were  compared  to 
forecasts  done  using  simple  exponential  smoothing  and  to  the 
actual  demand.  The  first  type  of  Box-Jenkins  model  built 
was  the  multivariate  (transfer  function)  model.  In  these 
models,  flying  hours  are  the  independent/input  variable  and 
demand  is  the  dependent/output  variable  for  forecasting. 

The  second  model  is  the  univariate  model  in  which  past  demand 
relationships  are  used  to  forecast  demand. 

The  three  models  forecast  one  quarter  of  demand.  The 
results  were  compared  to  the  actual  demand.  Results  showed 
low  correlation  between  flying  hours  and  demand  in  the  multi¬ 
variate  models.  Though  each  type  of  model  forecast  well, 
simple  exponential  smoothing  had  better  results.  In  the 
majority  of  forecasts,  the  three  models  overestimated  demand. 
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FORECASTING  RECOVERABLE  SPARES  USING 
BOX-JENKINS  TIME  SERIES  TECHNIQUES 


I .  Introduction 


General  Issue 

In  1982-83,  the  Air  Force  Chief  of  Staff  ordered  a 
comprehensive  study  of  the  aircraft  recoverable  spares 
forecasting  process  to  analyze  understated  forecasts 
totalling  millions  of  dollars.  An  aircraft  recoverable 
spare  is  a  part  that  is  normally  not  expended  in  use  and  can 
be  reused  after  recovery  and  repair  (8:571).  One  of  the  main 
areas  of  concern  for  Air  Staff  was  the  requirements  computa¬ 
tion  system  (7:58).  The  requirements  computation  system 
assumed  a  linear  causal  relationship  between  both  flying 
hours  and  demand  over  time.  This  linear  regression  technique 
known  as  Peacetime  Operating  Stock  Spares  Estimation  (POSSEM) 
underestimated  spares  demand  by  more  than  $1.3  billion  for 
1983  -  1985  (10:2).  Spares  budget  requirements  for  the 
future  were  determined  based  on  the  current  forecast  demand. 

Another  concern  was  the  increased  cost  of  weapon  system 
spares.  This  concern  has  magnified  the  search  for  improved 
forecasting  techniques  because  "reparable  spares  of  aircraft 
components  amounted  to  a  coirputed  budget  requirement  of 
$4,077  billion  is  FY  1985"  (14:v).  As  costs  increase,  recent 
cuts  in  the  Air  Force  budget  are  forcing  managers  to  more 
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accurately  predict  future  item  demand  for  recoverable  air¬ 
craft  spares.  Research  of  forecasting  techniques  is  ongoing. 

Countless  manhours  have  been  spent  by  the  Logistics 
Management  Institute,  the  Air  Force  Logistics  Command,  and 
the  rogistics  Mangement  Center  in  researching  and  testing 
current  forecasting  techniques  for  estimating  recoverable 
spares  demand.  Many  researchers  recommend  linear  regression 
and  forecasting  techniques  such  as  exponential  smoothing 
which  stress  linear  relationships.  Other  current  research 
has  questioned  the  assumption  that  the  relationship  between 
flying  hours  and  demand  is  linear.  Analyzing  the  question 
of  linearity  between  flying  hours  and  demand  is  a  key  step 
in  determining  which  forecasting  techniques  best  predict 
future  spares  requirements. 

Research  Question 

The  current  requirements  computation  system  understates 
recoverable  spares  demand  at  a  time  when  budget  constraints 
are  increasing.  Many  forecasting  techniques  in  use  base  the 
relationship  between  flying  hours  and  demand  as  linear  and 
do  not  consider  the  time  factor  in  both  demands  and  flying 
hours.  Emphasis  has  been  placed  on  finding  forecasting 
techniques  which  accurately  depict  recoverable  spares  demand. 
How  accurately  can  a  time  series  forecasting  model  using 
flying  hours  and  recoverable  spares  demand  predict  future 
recoverable  spares  requirements? 


Justification 


Many  forecasting  methods  compute  requirements  based  on 
the  linear  relationships  betveen  demand  and  flying  hours  with¬ 
out  considering  both  the  time  factor  and  non-linearity  involved 
in  demands  for  aircraft  recoverable  spares  and  flying  hours. 

In  the  thesis  "An  Alternative  Forecasting  Method  for  DRIVE," 
by  Captain  Alan  Closson,  he  discussed  an  example  of  a  model 
currently  in  use  by  the  Air  Force  which  assumes  recoverable 
spares  demand  varies  linearly  with  the  number  of  flying 
hours  (5:5).  The  model  is  the  Distribution  and  Repair 
in  Variable  Environments  (DRIVE)  developed  by  the  RAND 
Corporation.  DRIVE'S  purpose  is  to  ensure  repair  and  dis¬ 
tribution  of  recoverable  components  yields  the  best  support 
available  in  regenerating  spares  to  meet  demands  (6). 

The  DRIVE  algorithm  uses  a  four-quarter  moving  average 
to  determine  expected  demands  at  bases.  All  the  bases  are 
then  averaged  to  get  a  worldwide  average.  To  get  each 
bases  demand  rate,  the  worldwide  demand  rate  per  1,000  flying 
hours  is  multiplied  by  the  projected  number  of  flying  hours 
at  the  base.  Each  base's  demand  rate  is  based  on  the  assump¬ 
tion  demand  varies  linearly  with  the  number  of  flying 
hours  ( 6 ) . 

To  test  DRIVE'S  assumption,  the  Air  Force  Logistics 
Management  Center  ( LMC )  investigated  the  linear  relationship 
of  demand  and  flying  hours  since  forecasts  derived  from  world¬ 
wide  demands  are  based  on  linearity.  LMC  researchers  used 
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176  items  consumed  at  four  different  bases  with  varied  num¬ 
bers  of  aircraft.  The  results  of  their  study  indicated 
demand  does  not  vary  linearly  with  the  number  of  aircraft  at 
a  base  ( 1 ) . 

Two  other  examples  of  requirements  computation  methods 
currently  used  by  the  Air  Force  which  use  linear  causal  re¬ 
lationships  are  the  Peacetime  Operating  Stock  Spares 
Estimates  (POSSEM)  and  the  Air  Logistics  Early  Requirements 
Techniques  (ALERT)  which  were  studied  in  a  thesis  by 
Thomas  Lockette  (10).  Lockette's  research  also  indicated 
demand  and  flying  hours  are  not  just  linearly  related  (10:8). 
The  results  of  his  thesis  are  described  in  more  depth  in 
Chapter  II. 

A  technique  to  study  spares  forecasting  which  is  not 
based  on  linearity  is  time  series  analysis.  Time  series 
analysis  considers  factors  such  as  seasonality,  trends, 
cycles,  and  time.  Time  series  analysis  techniques,  such  as 
Box-Jenkins,  analyze  data  as  it  occurs  over  time  to  relate 
factors  like  seasonality,  trend,  and  cycles.  Seasonality  is 
a  condition  observed  when  a  systematic  pattern  occurs  in  time 
series  data.  An  example  is  high  demand  in  the  winter  and  low 
demand  in  the  summer  for  snow  shoes  (11:581).  Trend  is  a 
condition  observed  when  a  time  series  shows  consistent 
increase  or  decrease  over  time  (11:585).  A  cycle  is  a 
recurring  sequence  of  expansion  and  depression  periods. 
Examples  are  business  cycles  which  vary  in  length  and 
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magnitude  (11:567).  By  considering  non-linearity  over  time 
between  flying  hours  and  demand,  a  time  series  forecasting 
model  could  improve  future  requirements  computation  and  aid 
in  deciding  which  recoverables  to  budget  for  and  how  many  of 
each . 

Scope  and  Limitations 

Each  aircraft  weapon  system  is  composed  of  hundreds  of 
recoverable  spares.  Some  of  these  spares  can  be  used  on 
multiple  weapon  systems  while  others  are  unique  to  their 
aircraft.  Due  to  the  magnitude  of  aircraft  types  and  spares 
available,  the  scope  of  this  thesis  will  be  limited  to  ana¬ 
lyzing  C-135  aircraft  monthly  flying  hours  and  recoverable 
spares  demand  for  a  random  sample  of  spares  from  January  1978 
through  December  1985. 

A  large  number  of  time  series  analysis  and  forecasting 
software  are  available  with  which  to  analyze  data.  The  TIMES 
software  program,  which  was  developed  under  contract  for  the 
Air  Force,  and  FORECAST  MASTER  software  are  used  for  data 
analysis  and  testing  of  both  time  series  and  simple  exponen¬ 
tial  smoothing  forecasting  models.  The  TIMES  program  uses 
Fortran  subroutines  for  time  series  analysis  and  forecasting. 
The  program  is  available  on  the  Cyber  and  Harris  computers  at 
the  Air  Force  Institute  of  Technology  (AFIT).  It  was 
recently  adapted  for  use  on  personal  computers  at  AFIT. 

TIMES  does  Box-Jenkins  time  series  techniques  for  model 
building  (20).  FORECAST  MASTER  software  by  Scientific 
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Systems,  Electric  Power  Institute,  is  capable  of  doing 
limited  time  series  analysis,  regression,  and  simple  exponen¬ 
tial  smoothing  as  well  as  other  forecasting  techniques.  The 
times  series  and  exponential  smoothing  analyses  are  limited 
to  the  capabilities  of  each  of  these  software  packages. 


Investigative  Questions 

Answering  the  following  investigative  questions  will  aid 
in  the  analysis  and  development  of  the  time  series  and  simple 
exponential  smoothing  models  used  to  forecast  C-135  spares 
requirements . 

1.  What  type  of  statistical  relationship  exists  between 
C-135  flying  hours  and  demand  in  the  transfer  func¬ 
tion  models? 

a.  What  is  the  extent  of  correlation? 

b.  Are  there  linear  relationships  between  aircraft 
recoverable  demands  and  flying  hours? 

c.  How  are  past  C-135  spares  requirements  for  a 
specific  time  period  related  to  flying  hours  in 
one  or  more  prior  time  periods? 

d.  What  seasonality  or  trend  is  in  the  data? 

2.  Which  forecasting  techniques  best  model  the  data? 

a.  Do  the  transfer  function  models  simulate  the 
data  pattern? 

b.  What  time  series  parameters  were  found  in  the 
Autoregressive  Integrated  Moving  Average 

( ARIMA)  model? 

3.  How  accurately  do  the  developed  models  predict 
demand? 

a.  How  do  forecasts  from  the  transfer  function 
models  compare  with  actual  demand,  univariate 
models',  and  simple  exponential  smoothing 
models'  results? 

b.  What  modifications  are  needed? 

c.  Will  the  models  improve  spares  forecasting? 
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Methodology 

Three  types  of  models  are  built  to  compare  the  accuracy 
of  forecast  demand  to  actual  demand  for  ten  separate  C-135 
aircraft  recoverable  spares.  Box-Jenkins  time  series  tech¬ 
niques  are  used  to  develop  two  of  these  models  on  TIMES  soft¬ 
ware.  The  first  model  is  the  multivariate  or  transfer 
function  model  where  flying  hours  (independent  variable)  are 
used  to  predict  spares  demand.  The  second  Box-Jenkins  models 
developed  for  the  ten  spares  are  univariate  models  where  past 
demand,  indexed  over  time,  is  used  to  predict  future  demand. 
The  third  type  of  models  used  are  simple  exponential 
smoothing  models  developed  using  FORECAST  MASTER  software. 

Once  the  models  are  developed,  each  forecasts  one 
quarter  of  demand.  These  forecasts  are  then  compared  to  the 
actual  demand  for  the  quarter.  The  models'  forecasts  are 
compared  to  the  actual  values  by  calcualting  the  mean 
absolute  percentage  error. 

Summary 

Due  to  recent  budget  cuts  and  an  Air  Staff  directive, 
indepth  research  has  been  aimed  at  improving  the  forecast¬ 
ing  of  demand  for  aircraft  recoverable  spares.  Results  of 
research  indicate  that  the  relationship  between  flying  hours 
and  demand  may  not  be  linear.  This  thesis  will  investigate 
forecasting  C-135  recoverable  demand  by  using  Box-Jenkins 
time  series  techniques  on  TIMES  software.  In  Chapter  II,  a 
literature  review  is  given  on  Box-Jenkins  techniques  followed 
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by  a  reviev  of  results  found  in  three  Air  Force  studies  vhich 
implemented  Box-Jenkins  techniques  in  transfer  function  mo¬ 
dels.  Since  the  results  for  the  Box-Jenkins  models  will  be 


compared  to  those  found  using  exponential  smoothing,  a  reviev 
vill  also  be  given  on  simple  exponential  smoothing  and  some 
of  the  studies  researching  it.  To  conclude  Chapter  II,  a 
reviev  is  also  given  on  the  systems  American  Airlines  uses 
to  forecast  recoverable  aircraft  spares.  In  Chapter  III, 
the  methodology  for  each  step  of  the  Box-Jenkins  technique 
is  revieved  as  veil  as  the  data  conversion  steps  used  to 
prepare  the  spares  data.  Chapter  IV  is  an  analysis  and 
description  of  the  models  and  provides  a  comparison  of  their 
respective  forecasts.  Chapter  V  presents  the  conclusions  of 
the  research  and  suggests  areas  of  further  research. 
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I I .  Literature  Review 


Overview 

This  section  opens  with  a  brief  discussion  of  the 
Box-Jenkins  time  series  technique  for  developing  a  fore¬ 
casting  model.  This  is  followed  by  a  review  of  three  studies 
done  on  Air  Force  assets  which  used  Box-Jenkins  time  series 
techniques.  The  next  section  describes  the  simple  exponen¬ 
tial  smoothing  method  which  is  also  used  to  forecast 
spares  demand  so  a  comparison  can  be  made  between  the  results 
from  Box-Jenkins  techniques  and  simple  exponential  smoothing. 
The  final  section  describes  the  forecasting  methods  used  by 
American  Airlines. 

Box-Jenkins 

Box-Jenkins  is  a  time  series  forecasting  method  which 
uses  a  complex,  computer-based  procedure  to  produce  auto¬ 
regressive,  integrated  moving  average  (ARIMA)  models.  It  can 
adjust  for  seasonal  and  trend  factors  and  estimate  appropriate 
weighting  parameters  (9:123).  Box-Jenkins  consists  of  a  four- 
step  process  of  identification,  estimation,  diagnostic  check¬ 
ing  and  forecasting.  For  the  first  step  of  identification, 
George  Box  and  Gwilym  Jenkins  recommend  that  fifty  or  more 
observations  are  needed  before  their  correlogram  methods  pro¬ 
vide  a  reliable  guide  to  model  selection  for  monthly  seasonal 
data  such  as  spare  parts  demand  (2:64).  In  the  identification 
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step,  historical  data  are  analyzed  to  identify  a  model. 

Unknown  parameters  are  estimated  for  the  models  in  the  esti¬ 
mation  step.  During  diagnostic  checking,  the  models  are 
tested  for  validity  and,  if  valid,  they  are  used  in  the  fourth 
step  of  forecasting  where  future  values  are  likely  to  be 
observed  based  on  the  accuracy  of  the  developed  model  (2:64). 
Box-Jenkins  steps  are  used  to  develop  several  models  that 
explore  the  relationships  within  and  between  variables  such  as 
demand  and  flying  hours.  Once  the  initial  relationships  are 
found,  they  are  combined  into  a  final  transfer  function  model 
for  forecasting. 

Box-Jenkins  is  a  powerful  time  series  analysis  tool  that 
considers  relationships  between  variables  and  their  relation¬ 
ships  with  time.  Advantages  of  Box-Jenkins  are  that  it  allows 
a  wide  choice  of  weights  for  the  user  to  identify  much  more 
subtle  patterns  in  the  data  and  it  is  very  powerful  for  fore¬ 
casting  1  to  2  years  ahead  when  little  is  known  about  the 
underlying  pattern  of  the  data  (20:292).  Since  Box-Jenkins 
is  more  than  a  technique,  its'  philosophy  is  based  on  the 
principle  where  the  simpler  the  model,  the  better  it  is  so 
long  as  it  passes  diagnostic  checks  (11:98).  Disadvantages 
of  Box-Jenkins  are  that  it  is  complex  and  difficult  to  under¬ 
stand  and  it  is  expensive  in  its  use  of  computing 
time  (11:98).  Box-Jenkins  is  a  complex  technique  to  under¬ 
stand.  In  the  1970's,  it  was  usually  run  on  a  mainframe 
computer  which  was  time  consuming  and  expensive.  Today, 
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software  exists  that  runs  complex  methods  such  as  Box-Jenkins 
on  personal  computers  allowing  a  hands  on  approach  to 
computer  analysis  with  less  expense. 

Research  Studies  Using  Box-Jenkins  Techniques 

Three  studies  using  Box-Jenkins  time  series  techniques 
were  accomplished  recently  by  Air  Force  personnel.  The  first 
was  done  in  1982  by  Singpurwalla  and  Talbott  (15:552).  They 
used  time  series  analysis  to  investigate  the  interrelation¬ 
ships  between  alert  availability  and  flying  hours  for  Air 
Force  C-141  aircraft.  Singpurwalla  and  Talbott  used  the  Box- 
Jenkins  method  to  develop  a  transfer  function  model.  Using 
the  final  transfer  function  model,  their  results  indicated 
flying  hours  do  have  an  effect  on  the  alert  availability 
specified  by  the  functional  form  of  their  developed  final 
transfer  function  model  (15:582).  Singpurwalla  and  Talbott 
concluded  their  final  equation  supported  the  expression 
".  .  .  the  more  you  fly,  the  less  you  fail"  (15:582). 

In  1983,  Larry  Taylor  did  a  masters  thesis  on  using 
time  series  analysis  to  develop  forecasts  for  recoverable 
spares.  He  used  five  F-16  line  replacement  units  (LRU) 
that  impaired  mission  capability  (MICAP)  for  the  week  of 
12  June  1983  (16:8).  He  wanted  to  demonstrate  that  the  num¬ 
ber  of  demands  per  flying  hour  is  not  a  linear  ratio. 

Taylor  chose  time  series  analysis  because  it  does  not  rely 
on  cumulative  averages,  but  uses  actual  data  as  it  occurs 
through  time  to  make  forecasts  (16:8).  Taylor  not  only  noted 
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that  the  "linearity"  assumption  caused  support  problems,  he 
also  noted  that  demands  per  flying  hour  data  are  collected 
over  time  so  time  is  also  a  factor.  Since  time  is  a  factor, 
there  are  underlying  interrelationships  between  demand  and 
flying  hours  (16:28).  Time  series  analysis  is  a  viable 
approach  to  help  determine  if  those  relationships  are  linear 
or  not . 

Taylor  used  the  TIMES  software  package  for  analysis.  In 
his  findings,  Taylor  noted  a  seven  period  time  lag  for  demand 
before  flying  hours  had  an  effect  on  it.  It  was  also  shown 
both  flying  hours  and  demand  were  nonstationary  and  had  auto¬ 
regressive  and  moving  average  components  indicating  possible 
non-linearity  (16:69-70).  Taylor  recommended  further  study 
since  he  used  less  than  50  observations  and  recommended  using 
at  least  50  and  preferably  100  observations  to  realize  the 
full  benefit  of  Box-Jenkins  time  series  analysis  (16:72). 

In  a  third  study  using  Box-Jenkins,  Thomas  Lockette 
did  a  masters  thesis  in  1984  to  accomplish  a  time  series 
analysis  of  C-130  recoverable  spares  requirements  to 
identify  the  relationship  between  C-130  flying  hours  and 
spares  requirements.  The  study  was  accomplished,  not  to 
develop  a  forecasting  model,  but  to  demonstrate  a  methodology 
to  be  used  to  identify  relationships  for  various  logistics 
times  series  for  other  weapon  systems  (10: ii). 

Lockette  examined  the  relationships  between  flying  hours 
and  expenditures.  His  research  was  done  to  investigate  how 
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Box-Jenkins  time  series  techniques  could  be  used  to  develop 
models  which  would  forecast  Program  Objective  Memorandums 
(POMs).  He  then  compared  his  formula  to  the  Peacetime  Oper¬ 
ating  Stock  Spares  Estimation  System  (POSSEM)  linear  regres¬ 
sion  technique  (10:4).  Lockette  attempted  to  demonstrate 
that  current  linear  regression  techniques  used  to  forecast 
obligations  could  be  improved  by  using  Box-Jenkins  techniques 
to  recognize  time  factors  inherent  in  the  data  (10:75).  He 
obtained  conflicting  results  and  mentions  the  need  for  fur¬ 
ther  study  using  demand  versus  flying  hours  instead  of  ob¬ 
ligations  and  flying  hours.  He  recommended  using  well  over 
50  data  points  since  he  could  not  confirm  if  his  results, 
indicating  the  presence  of  some  non-linearity,  might  have 
resulted  from  using  a  small  number  of  data  points  (10:76). 

Over  each  of  these  studies,  the  authors  mentioned  their 
results  using  Box-Jenkins  time  series  techniques  were  limited 
due  to  the  data  constraints.  This  research  attempts  to  build 
on  these  three  previous  studies  by  developing  both  the  multi¬ 
variate  and  univariate  Box-Jenkins  models  to  forecast 
aircraft  demand  and  then  comparing  them  to  a  currently  used 
method  of  forecasting,  simple  exponential  smoothing.  Since 
the  Box-Jenkins  models’  forecasts  are  compared  to  simple 
exponential  smoothing,  the  next  section  reviews  simple 
exponential  smoothing  and  research  results  found  using  it. 
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Exponential  Smoothing 

Exponential  smoothing  is  a  type  of  moving  average  where 
past  data  points  are  not  given  equal  weight.  Recent  data  is 
weighted  more  heavily  than  older  data.  The  value  of  a 
smoothing  constant  determines  the  level  of  weight  assigned  to 
most  recent  data.  The  smoothing  constant  takes  on  a  value 
between  one  and  zero,  where  small  values  smooth  recent  trends 
and  large  values  emphasize  recent  demand  conditions  (18:53). 

The  formula  for  single  exponential  smoothing  is: 

S  =  aX  +  (l-a)S  (1) 

t  +  1  t 


where : 

X  =  most  recent  actual  value 

S  =  the  latest  forecast 
t 

3  =  forecast  for  the  next  period 

t  +  1 

a  =  the  smoothing  constant  (20:79). 

Simple  exponential  smoothing  is  an  extremely  popular 
technique  among  forecasting  practioners  since  the  logic  of 
the  whole  exponential  smoothing  process  is  based  on  the 
following  premise:  "If  the  forecast  for  a  particular  period 
was  too  high,  reduce  it  for  the  next  period;  if  it  was  too 
low,  raise  it!"  (11:107).  Advantages  of  exponential  smooth¬ 
ing  are  that  it  is  easy  to  computerize  for  a  large  number  of 
items  at  low  cost,  easy  to  monitor,  and  easy  to  under¬ 
stand  (11:98).  The  major  advantage  of  exponential 


smoothing  is  that  it  does  not  require  a  long  history  of  data 
to  get  started.  It  captures  the  effect  of  all  past  data 
on  the  most  previous  forecast  (18:53).  Disadvantages  of 
exponential  smoothing  mentioned  by  Makridakis  and  Wheel¬ 
wright  are;  it  lags  shifts  in  demands,  it  can  only  be 
applied  to  short-run  forecasts,  it  may  overreact  to  random¬ 
ness,  and  technically,  problems  arise  in  choosing  smooth¬ 
ing  constant  values  (11:110).  Several  studies  have  been 
done  using  exponential  smoothing  to  analyze  aircraft  spares 
demand  data  for  the  Air  Force. 

In  a  study  of  Air  Force  spare  parts  done  in  1984, 

Craig  Sherbrooke  analyzed  worldwide  base  level  demands  for 
1,027  recoverable  spares  over  48  months.  He  found  that 
demand  in  consecutive  months  was  not  independent  since  many 
spares  had  mean  demands  that  were  not  constant.  He  reported 
a  model  would  be  needed  that  would  not  assume  demand  through 
time  is  independent.  Sherbrooke  concluded  exponential 
smoothing  was  a  better  predictor  of  mean  demand  since  mean 
demand  rates  change  over  time  and  exponential  smoothing  will 
be  more  responsive  to  change  (13:24). 

In  1987,  Sherbrooke  did  another  study  on  aircraft  recov¬ 
erable  spares.  He  compared  the  Air  Force  Recoverable  Con¬ 
sumption  Item  Requirements  (D041)  system  moving  average 
rates  to  those  found  using  exponential  smoothing  and  Bayes 
predictions.  Demand  data  were  used  from  560  C-5  aircraft 
first  indenture  non-engine  items  over  16  quarters.  This 


data  was  accumulated  from  the  D041  system.  Demand  data  for 
the  A-10  airframe  and  F-16  engine/airframe  were  also  used 
in  further  testing.  Two  important  results  Sherbrooke  found 
were:  (1)  demand  in  any  period  is  more  highly  correlated 
with  recent  demand  than  with  earlier  demand,  (2)  exponential 
smoothing  techniques  are  superior  to  today's  D041  moving 
average  rates  for  estimating  demand  rates  (14:11).  In  both 
his  studies,  Sherbrooke  recommends  a  smoothing  constant  value 
of  0.4  because  it  reduced  mean  squared  error  and  mean 
absolute  error  for  the  quarterly  data  (13:23). 

A  final  step  in  searching  the  literature  was  to  look  for 
information  on  how  United  States  commercial  airlines  forecast 
their  recoverable  spares.  What  follows  is  a  brief  summary 
of  how  American  Airlines  forecasts  rotable  requirements. 
Rotable  is  the  American  Airlines  term  for  recoverable 
assets  (17:1). 

American  Airlines 

Air  Force  C-135  aircraft  fly  missions  ranging  from  cargo 
and  passenger  transport  to  aerial  refueling.  They  are  physi¬ 
cally  structured  similar  to  commercial  aircraft  and  also  have 
similar  maneuverability  and  flight  capabilities. 

Due  to  the  similarities  of  Air  Force  C-135's  to  com¬ 
mercial  aircraft,  requests  were  made  to  American  Airlines  for 
information  on  how  they  forecast  recoverable  spares. 

Mark  Tedone,  American  Airlines  Decision  Technologies  (AADT), 
has  written  an  article  for  Interfaces  titled  "Reparable 
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Part  Management"  which  is  scheduled  to  be  published  some  time 
in  1989.  His  article  describes  the  new  Rotables  Allocation 
Planning  System  (RAPS)  recently  developed  by  AADT  H.  .  .to 
provide  forecasts  of  rotable  parts  demand,  recommend  least-cost 
allocations,  and  calculate  the  availability  level  associated 
with  the  optimal  allocation  of  each  part”  (17:ii).  The  pre¬ 
decessor  to  RAPS  is  the  Rotable  Forecasting  and  Availability 
Control  System  (ROFACS).  It  was  developed  in  the  late  1960's 
and  implemented  in  the  mid  1970's.  It  was  originally  deve¬ 
loped  for  a  rigid  mainframe  computer.  ROFACS  generated  demand 
forecasts  using  only  time  series  methodology.  For  its  time, 
ROFACS  remained  a  valuable  indicator  of  appropriate  alloca¬ 
tion  levels.  Since  ROFACS  used  only  time  series  methodology 
on  a  rigid  mainframe  computer,  it  was  sometimes  slow  in  re¬ 
sponding  to  changes  in  aircraft  utilization  or  major  fleet 
expansions.  Decision  Technologies  examined  ROFACS  in  detail 
to  uncover  its  deficiencies  and  started  development  of 
RAPS  (17:3-4).  RAPS  uses  linear  regression  to  establish  a 
relationship  between  parts  removal  and  various  functions  of 
monthly  flying  hours  to  determine  the  absolute  volume  of 
demand  (17:6-7).  American  Airlines  also  determines  the 
patterns  of  demand  through  time  using  RAPS,  which  has  a  major 
influence  on  the  total  cost  of  an  allocation  (17:8).  Overall, 
American  Airlines  has  made  multi-million  dollar  improvements 
in  their  inventory  position  by  taking  their  older  time  series 
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methodology  for  a  mainframe,  finding  its  deficiencies,  and 
using  them  to  develop  a  new  RAPS  system  which  has  exceeded 
expectations  (17:15-16). 

American  Airlines  flies  a  city  hub  system  for  aircraft 
and  has  set  repair  sites  and  flight  routes.  Their  aircraft 
spend  few  hours  sitting  idle  since  a  grounded  aircraft  lowers 
customer  service.  American  Airlines  used  the  ROPACS  time 
series  approach  to  forecast  spares  when  mainframe  computers 
were  all  that  was  available.  They  have  recently  changed 
to  RAPS  and  now  use  linear  regression  techniques  on 
real-time  personal  computer  systems  throughout  their  loca¬ 
tions  to  forecast  spares.  RAPS  is  also  used  for  recommend¬ 
ing  least-cost  allocations  of  spares  and  calculating 
availability  levels  for  optimum  spares  allocation.  The 
RAPS  implementation  provided  a  multi-million  dollar  bene¬ 
fit  to  American  Airlines  and  although  they  currently  use 
linear  regression  without  time  series  methods,  it  was  a 
time  series  methodology  which  broke  ground  for  the  RAPS. 

Summary 

Chapter  II  provided  an  introduction  to  Box-Jenkins 
techniques  and  the  four-step  process  of  identification, 
estimation,  diagnostic  checking,  and  forecasting.  A  review 
was  given  of  previous  Box-Jenkins  studies  done  on  Air  Force 
assets.  After  Box-Jenkins,  simple  exponential  smoothing  was 
reviewed  as  well  as  studies  done  by  Craig  Sherbrooke  who 
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used  simple  exponential  smoothing  on  Air  Force  assets. 

The  chapter  concluded  with  a  discussion  of  how  American 
Airlines  forecasts  recoverable  spares. 

In  the  Chapter  III  methodology,  each  of  the  Box-Jenkins 
analysis  steps  used  to  answer  the  investigative  questions  and 
lead  to  the  development  of  a  model  will  be  discussed  as  well 
as  the  steps  taken  to  obtain  and  prepare  the  flying  hour  and 
demand  data  for  use  with  both  Box-Jenkins  and  simple  exponen¬ 
tial  smoothing  models. 
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III.  Methodology 


Overview 

This  chapter  opens  with  a  description  of  the  spares  and 
flying  hours  data  collection  and  preparation  for  Box-Jenkins. 
The  next  section  describes  the  Box-Jenkins  steps  taken  in 
using  historical  data  to  develop  a  model  for  forecasting. 

The  steps  of  model  identification,  parameter  estimation,  and 
forecasting  are  the  principal  stages  in  the  TIMES  computer 
program  that  will  be  used.  To  avoid  duplication  of  explana¬ 
tion,  a  general  description  will  be  given  of  each  step  since 
the  detailed  model  formulation,  plots,  and  tables  will  be 
shown  in  detail  in  Chapter  IV.  The  final  section  briefly 
describes  the  software  used  to  model  the  spares  data  using 
exponential  smoothing. 

Data  Availability  and  Preparation 

Flying  Hours .  Monthly  C-135  flying  hour  data  were 
obtained  from  the  Air  Force  Inspection  and  Safety  Center  at 
Norton  AFB,  CA,  for  the  period  of  January  1978  -  March  1986. 
Ninety-six  months  were  used  for  model  building.  Three 
months  of  data  were  withheld  as  a  basis  to  compare  the  fore¬ 
casts  with  the  actual  observations.  The  data  are  attached 
in  Appendix  A. 
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Recoverable  Spares .  The  recoverable  spares  demand 
data  for  C-135  aircraft  was  obtained  from  the  Analysis  and 
Information  Management  Division  ( MMMA )  at  Headquarters, 

Air  Force  Logistics  Command.  The  quarterly  data  were  down¬ 
loaded  from  a  tape  containing  D041  data  from  January  1978  - 
March  1986.  Ten  C-135  spares  were  randomly  chosen  from  a 
listing  of  200  recoverable  spares'  data  downloaded  by 
HQ  AFLC/MMMA.  Table  1  lists  the  stock  numbers  and  their 
description  as  they  appeared  on  the  D041.  Throughout  the 
thesis,  the  stock  numbers  are  identified  by  the  abbrevia¬ 
tions  in  the  third  column  of  Table  1.  The  quarterly 
demand  data  are  attached  in  Appendix  A.  Of  the  thirty-three 
quarters  of  data  for  each  spare,  seven  quarters  were  missing. 
As  with  the  flying  hours,  the  quarter  ending  March  1986  was 
withheld  as  a  basis  for  forecast  comparisons.  Asterisks 
appear  in  Appendix  A  next  to  the  dates  of  the  missing 
quarters  of  demand  data. 


Stock  Number 

1560  00  075  1401 
1560  00  056  6631 
1560  00  011  0674 
1560  00  015  2153 
1560  00  055  5584 
1560  00  055  5583 
1560  00  055  5585 
1560  00  056  6632 
1560  00  056  6633 
1560  00  068  2535 


Table  1. 

-135  Spares'  Stock 

Description 

FLFAIRING 
FLPANEL  ASSY 
FLTAB  ASSY 
FLFRAME 
FLTAILPIPE 
FLCOWL 
FLSLEEVE 
FLPANEL  ASSY 
FLPANEL  ASSY 
FLFORK  ASSY 


Numbers 

Abbreviations 

FAIR01 

ASSY31 

TAB7  4 

FRAME53 

TAIL84 

COWL 8 3 

SLVE85 

ASSY32 

ASSY33 

FORK35 


21 


To  find  values  for  the  missing  quarterly  data,  an  Impu¬ 
tation-based  procedure  was  used  known  as  mean  imputation 
where  means  from  sets  of  values  are  used  to  fill  in  missing 
data  (12:6).  To  obtain  the  sets  of  values  for  each  spares' 
seven  missing  quarters,  a  modified  Delphi  technique  was  used. 
Three  supply  experts  separately  estimated  the  demand  for  each 
missing  quarter  by  3tock  number  by  analyzing  the  pattern  of 
available  quarterly  data.  The  mean  of  the  experts'  deter¬ 
minations  was  used  as  the  demand  for  each  of  the  missing 
quarters  of  data. 

For  data  analysis,  32  of  the  33  quarters  were  used  for 
model  building.  The  remaining  quarter  (Jan  -  Mar  1986), 
was  withheld  for  comparison  with  forecasts  from  the  models. 

To  match  the  96  months  of  flying  hours  data,  the  32 
quarters  of  data  were  converted  to  96  months  by  dividing 
each  quarter's  spares  demand  by  3 .  In  a  Sherbrooke  study, 
he  briefly  mentioned  converting  quarterly  data  to  monthly 
is  feasible  since  it  induces  an  insignificant  amount  of 
error.  Overall,  the  amount  of  error  induced  depends  totally 
on  the  distribution  of  the  data  and  the  magnitude  of  the 
values  (13:10) . 

Box-Jenk ins  Model  Development 

Box-Jenkins  time  series  analysis  uses  a  four-step  process 
of  model  identification,  parameter  estimation,  diagnostic 
checking,  and  forecasting.  These  four  steps  are  divided  into 
three  stages.  Initially,  a  general  class  of  forecasting 
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models  Is  postulated.  In  stage  1,  a  specif Lc  model  is  tenta¬ 
tively  picked  as  the  best  forecasting  method  for  that  set  of 
data.  In  stage  2,  the  model  is  fitted  to  the  historical  data 
and  tested.  If  it  is  not  adequate,  stage  1  is  repeated.  Once 
the  model  passes  the  checks,  it  is  used  to  develop  a  forecast 
in  stage  3  (19:173). 

Model  Identification 

In  the  model  identification  stage,  the  historical  data 
for  flying  hours  and  demand  are  separately  analyzed  to  iden¬ 
tify  the  relationships  that  exist  within  each  time  series. 
Identifying  possible  models  involves  the  concepts  of  station- 
arity,  correlation,  autocorrelation,  and  partial  autocorrela¬ 
tion.  First,  to  develop  each  individual  univariate  time 
series  model,  each  set  of  data  must  be  stationary.  A  data 
set  is  considered  to  be  stationary  if  its  values 
(observations)  fluctuate  around  a  constant  mean  with  a  con¬ 
stant  variance.  Stationarity  is  checked  by  observing  plots 
of  the  original  data  and  residual  plots  output  by  the  TIMES 
computer  program.  If  a  model  is  non-stat ionary,  two  commonly 
practiced  techniques  for  making  it  stationary  are  to  trans¬ 
form  the  original  data  into  the  natural  logarithm  of  the 
original  or  to  difference  the  data.  First  differencing 
of  data  is  employed  when  just  the  mean  is  changing.  Second 
differencing  is  usually  done  when  both  the  mean  and  slope 
are  changing.  First  differencing  is  accomplished  by 
subtracting  successive  values  from  one  another  and  using 
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these  values  as  the  nev  data  set.  The  first  differenced 
series  has  one  less  point  than  the  original  data.  Second 
differencing  is  accomplished  using  the  first  differenced  data 
set,  differencing  it  again,  and  then  using  that  series.  A 
second  differenced  data  set  has  two  less  data  points  than  the 
original  data  (4).  Once  the  data  are  stationary,  development 
of  the  models  takes  place. 

In  developing  each  model,  correlation,  autocorrelation, 
and  partial  correlation  need  to  be  known.  Correlation  is  the 
association  or  mutual  correspondence  between  two  variables 
and  describes  what  happens  to  one  of  the  variables  if  there  is 
a  change  in  the  other  (19:174).  Autocorrelation  describes  the 
association  or  mutual  dependence  of  the  same  variable  at 
different  time  periods.  Autocorrelation  is  critical  since 
" .  .  .  autocorrelation  among  successive  values  of  the  data 
is  a  key  tool  in  identifying  the  basic  pattern  (or  rather  the 
model  corresponding  to  it)  that  describes  the  data"  (19:174). 
Partial  autocorrelat ions  are  analogous  to  autocorrelat ions 
because  they  indicate  relationships  of  the  time  series  val¬ 
ues  to  various  time  lagged  values  of  the  same  series.  They 
differ  from  autocorrelation  since  they  are  computed  for  each 
time  lag  after  the  effect  of  all  other  time  lags  on  the  given 
one,  and  on  the  original  data  are  removed  (16:30). 

The  autocorrelation  is  used  to  help  identify  the  time 
series  data  pattern,  test  for  stat ionar i ty,  and  identify  a 
model  for  the  time  series.  The  autocorrelation  coefficient 
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value  which  measures  the  linear  relationship  between  two  time 
series  observations  is  always  between  -1  and  +1.  If  the  value 
is  near  0  (zero),  the  two  observations  are  linearly  indepen¬ 
dent.  If  the  value  is  near  -1  or  +1,  the  two  observations 
are  more  linearly  dependent  (10:17).  Once  each  model  is 
stationary  and  all  the  autocorrelation  coefficients  are 
known,  the  autocorrelation  coefficients  are  converted  to  an 
autocorrelation  function  ( ACF ) .  Partial  autocorrelation 
functions  (PACF)  are  also  determined.  By  looking  at  the 
patterns  of  the  ACF  and  PACF,  a  tentative  model  is  developed 
which  could  indicate  the  presence  of  autoregressive  (AR)  or 
moving  average  (MA)  parameters  for  an  Autoregressive 
Integrated  Moving  Average  (ARIMA)  model  (19:178). 

Once  the  models  for  both  flying  hours  and  demand  are 
determined  to  be  either  AR  or  MA,  or  a  combination  of  both, 
they  can  be  written,  for  example,  as  either  an  ARIMA  (p,d,q) 
or  an  ARMA  (p,q)  in  Box-Jenkins  notation.  The  number  of 
significant  spikes  in  the  autocorerelat ion  and  partial  auto¬ 
correlation  functions  are  denoted  by  the  letters  p  and  q. 

The  d  represents  the  level  of  differencing  where  zero  means 
no  differencing  was  used.  Now  that  the  models  are  identi¬ 
fied,  the  next  step  is  parameter  estimation. 

Parameter  Estimation 

After  identifying  the  tentative  model,  the  parameter 
values  for  AR  or  MA  or  both  are  estimated  from  the  data. 

The  TIMES  program  is  run  on  an  AFIT  persona)  computer  to 
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obtain  the  parameter  estimates  which  minimize  the  mean  squared 
errors  (21:12).  Once  the  parameters  are  estimated  for  the 
tentative  model,  diagnostic  checks  are  performed  to  test  the 
model  adequacy. 

Diagnostic  Check ina 

"A  model  is  considered  adequate  if  the  residual  differ¬ 
ences  between  the  individual  time  series  observations  and  the 
forecast  using  the  tentative  model  are  white  noise"  (10:21). 
White  noise  means  the  residuals  show  no  linear  relationships 
between  themselves  so  they  cannot  be  forecast  by  an  improved 
model.  Box  and  Jenkins  emphasize  several  checks  that 
should  be  done  to  test  model  adequacy.  The  first  check 
strongly  emphasized  by  Box  and  Jenkins  is  to  visually  inspect 
the  plot  of  the  residuals  themselves  (3:289).  In  these  plots, 
the  mean  should  be  near  zero  and  should  not  be  changing  over 
time.  The  residual  variance  should  not  change  ever  time 
either . 

The  second  check  recommended  is  an  autocorrelation  check 
on  the  residual  ACF  and  PACF  to  make  sure  there  are  no  signi¬ 
ficant  spikes  which  could  indicate  more  AR  or  MA  parameters  are 
needed  (3:290).  A  third  check  for  model  adequacy  is  the  Port¬ 
manteau  lack  of  fit  test.  This  test  involves  solving  for  a 
variable  Q  and  comparing  it  to  Chi-Square  values  at  selected 
confidence  intervals  with  the  same  degrees  of  freedom.  If 
the  Portmanteau  Q-value  is  less  than  the  Chi-Square  value  at 
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the  same  degrees  of  freedom,  the  model  surpasses  the  tested 
confidence  level  (3:292).  The  fourth  test  to  check  is  model 
Inadequacy  arising  from  changes  in  parameter  values.  The 
parameter  values  should  remain  constant  over  time  (3:293). 

Four  other  diagnostic  checks  used  to  check  a  model  were 
discussed  in  LOGM  630,  Forecasting  Management,  by 
Lieutenant  Colonel  Bruce  P.  Christensen  (4).  The  four  other 
checks  are: 

1.  Cumulative  Periodogram.  This  periodogram  is  a  Fourier 
transform  of  error  covariance  and  is  an  output  from 
the  TIMES  computer  program.  Ideally,  the  periodogram 
will  show  a  diagonal  line  through  a  rectangle  which 
indicates  pure  white  noise  and  ensures  additional 
seasonality  is  not  in  the  residuals  (4). 

2.  Histogram  of  the  Residuals.  The  residual  histogram 
indicates  if  residuals  are  distributed  normal  or  not. 
If  they  are  not  normal,  the  model  can  still  be  good 
but  care  must  be  taken  in  giving  confidence  levels. 

The  histogram  is  also  an  output  of  TIMES  (4). 

3.  Variance  of  Residuals.  The  smaller  the  value  of 
the  variance  of  residuals,  the  better  the  model. 

The  variance  of  residuals  is  used  to  compare  alter¬ 
native  models.  The  better  model  is  usually  the  one 
with  the  lowest  residual  mean  square.  This  value  is 
calculated  and  also  output  by  TIMES  (4). 

4.  Power  Spectrum.  The  power  spectrum  shows  whether 
or  not  additional  seasonality  exists  in  the 
residuals.  The  power  spectrum  is  a  Fourier  trans¬ 
formation  of  autocovariance.  Ideally,  it  is  a 
smooth  horizontal  line  across  a  rectangle  (4). 


Once  each  model  for  flying  hours  and  demand  pass  the 
diagnostic  checks,  they  are  considered  prewhitened  and  are 
ready  to  be  combined  into  a  two  variable  (multivariate) 
model  called  the  transfer  function  model.  The  transfer 
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function  model  must  pass  through  the  same  steps  of  identifi¬ 
cation,  parameter  estimation,  and  diagnostic  checking  as  the 
univariate  models  already  discussed. 

Transfer  Function  Model 

As  univariate  models,  the  time  series  of  flying  hours 
and  demand  were  separately  analyzed  to  identify  the  "within" 
relationships  in  each  time  series.  In  the  transfer  func¬ 
tion  model  phase,  both  the  flying  hours  and  demand  models 
are  prewhitened  and  the  residuals  from  the  two  series  are 
analyzed  to  identify  the  "between"  relationships  of  the  two 
models.  The  two  variable  model  is  a  transfer  function  model. 

The  transfer  function  model  now  consists  of  residuals 
from  the  prewhitened  flying  hours  and  demand  models.  The 
residuals  are  analyzed  to  identify  the  relationships  which 
exist  between  the  two  time  series.  Flying  hours  is  the 
independent  univariate  and  demand  is  the  dependent  univariate. 
Noise  parameters  may  also  be  added  to  the  transfer  function 
model  to  identify  noise  which  is  not  white  noise  (totally 
random)  left  by  the  two  univariate  models.  Together  they 
form  the  final  transfer  model  used  to  forecast  demand. 

The  general  form  of  the  transfer  function  model  is: 

-1 

Y  =  s  ( B ) w ( B ) X  +  N  (2) 

t  t-1  t 
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(3) 


-1 

N  =  0  (B)0(B)a 

t  t 

where : 

a  =  white  noise 
t 

s,0,and  w  =  parameter  values 

X  =  prewhitened  independent/input  series 
t  (flying  hours) 

N  =  noise 
t 

Y  =  prewhitened  dependent/output  series 
t  (spares  demand)  (3:362). 

Exponential  Smoothing 

The  ten  separate  spares  models  were  built  using  the 
simple  exponential  smoothing  program  on  FORECAST  MASTER 
personal  computer  software  by  Scientific  Systems,  Electric 
Power  Research  Institute,  1986.  Each  model  was  built  with 
an  a  =  0.4  as  recommended  by  Sherbrooke.  Each  model  was 
then  used  to  forecast  its  respective  demand  for  the  period 
January  -  March  1986. 

Compar ison  of  Mode  1  Forecasts 

Once  the  Box-Jenkins  transfer  function,  Box-Jenkins 
univariate,  and  simple  exponential  smoothing  models  are 
developed  and  tested,  each  model  is  used  to  forecast  demand 
for  the  Jan  -  Mar  1986  quarter  for  each  of  the  spares.  To 
compare  the  forecasting  accuracy  of  the  three  techniques. 
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the  mean  absolute  percentage  error  is  calculated  for  each 
forecast.  The  mean  absolute  percentage  error  is  calculated 
by  subtracting  the  forecast  demand  value  from  the  actual 
demand  value  and  then  the  absolute  value  is  taken.  The 
absolute  value  is  then  divided  by  the  actual  value  and 
multiplied  Dy  luO  to  make  it  a  percentage.  The  mean  absolute 
percentage  error  was  chosen  since  it  is  used  to  evaluate  and 
compare  a  single  observation  forecast  by  the  three 
different  models  (4).  In  the  tables  shown  in  Chapter  IV, 
the  mean  absolute  percentage  error  is  denoted  as  the  MAPE. 

Investigative  Questions  Research 

As  the  models  are  built  and  used  to  forecast,  the  re¬ 
sults  are  used  to  answer  the  investigative  questions  asked  in 
Chapter  I.  Building  and  analyzing  the  transfer  function 
models  leads  to  an  answer  for  investigative  question  1 
which  is: 

1.  What  type  of  statistical  relationship  exists 
between  C-135  flying  hours  and  demand  in  the 
transfer  function  models? 

Investigative  questions  2  and  3  are  answered  once  the  models 
are  built  and  the  results  of  the  models'  forecasts  are 
compared.  Questions  2  and  3  are: 

2.  Which  forecasting  techniques  best  model  the 
data? 

3.  How  accurately  do  the  developed  models 
predict  demand? 
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Summary 

The  Box-Jenkins  four-step  process  of  model  identifica¬ 
tion,  parameter  estimation,  diagnostic  checking,  and  fore¬ 
casting  was  used  to  develop  multivariate  models  that  related 
the  relationships  within  and  between  flying  hours  and  demand 
to  forecast  future  spares  demand.  The  four-step  process  was 
also  used  to  build  univariate  demand  models  which  forecast 
demand  based  on  the  pattern  and  relationships  in  the  past 
demand.  The  simple  exponential  smoothing  models  based  future 
spares  demand  on  recent  past  demand  in  determining  the  fore¬ 
cast.  Both  the  univariate  models  and  simple  exponential 
smoothing  models  were  developed  to  compare  forecasting 
accuracy  with  the  transfer  function  models.  In  Chapter  IV, 
a  step-by-step  description  of  the  Box-Jenkin  four-step 
process  is  given  with  TIMES  output  shown  to  aid  in  analyzing 
the  models.  The  results  of  both  the  Box-Jenkins  and  simple 
exponential  smoothing  forecasts  are  compared  to  the  actual 
spares  demand. 
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IV.  ANALYSIS  AND  RESULTS 


Introduction 

This  chapter  describes  and  analyzes  the  steps  used  in 
developing  each  model  and  provides  a  comparison  of  the  models 
forecasting  results  for  each  of  the  Box-Jenkins  multivariate 
(transfer  function  models),  Box-Jenkins  univariate,  and  simple 
exponential  smoothing  spares  forecasting  models  for  data 
spanning  1978  -  1985.  The  chapter  opens  with  a  discussion  of 
each  Box-Jenkins  step  u.ied  to  form  the  univariate  flying 
hours  and  spares  time  series  models  which  were  then  used  to 
form  the  transfer  function  models  for  forecasting.  The 
section  concludes  with  an  analysis  of  the  multivariate 
forecasting  results. 

Second,  a  detailed  example  is  given  of  the  steps  taken 
in  buiding  one  of  the  univariate  Box-Jenkins  models  where 
past  demand  data  are  used  to  predict  future  demands.  In  the 
multivariate  models,  the  relationships  between  and  within 
flying  hours  and  demands  were  used  to  build  a  model  to  fore¬ 
cast.  A  discussion  of  the  univariace  models'  forecasting 
results  is  given  at  the  end  of  this  section. 

Third,  the  results  of  forecasting  demand  using  simple 
exponential  smoothing  is  shown  for  the  ten  spares.  The 
chapter  concludes  with  a  comparison  of  the  forecasting 
results,  a  discussion  of  the  answers  found  for  the  inves¬ 
tigative  questions  asked  in  Chapter  I,  and  a  summary. 
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Transfer  Function  Model  Development 

To  form  the  multivariate  time  series  models  for  each  of 
the  ten  spares,  the  flying  hours  data  and  spares  data  were 
used  to  form  prevhitened  univariate  models.  They  were  then 
combined  to  form  the  ten  transfer  function  models  used  to 
forecast  spares  demand.  A  description  of  the  steps  used  in 
developing  the  univariate  model  for  the  independent  variable 
(flying  hours)  is  described  below. 

Flying  Hours  Identification.  Before  the  transfer 
function  model  can  be  formed,  each  individual  series  must  be 
prewhitened.  The  identification  step  is  the  first  step  toward 
prewhitening  the  flying  hours  data.  The  flying  hours  data  for 
96  months  (Appendix  A)  were  analyzed  using  TIMES  and  FORECAST 
MASTER  software.  Prior  to  analysis,  the  logarithmic  function 
of  the  flying  hours  data  was  taken  because  the  flying  hours 
consisted  of  very  large  monthly  numbers  compared  to  monthly 
demand  for  the  ten  spares.  Taking  the  logarithmic  function 
of  the  flying  hours  reduces  the  magnitude  of  the  residuals  so 
impulse  response  weights  can  be  determined  during  develop¬ 
ment  of  the  transfer  function  model. 

To  start,  FORECAST  MASTER  was  used  to  obtain  a  plot  of 
the  transformed  data  (Figure  1)  and  TIMES  was  used  to 
perform  an  initial  identification  of  data  stationar ity . 

The  ACF  did  not  drop  rapidly  to  zero  and  the  TIMES  run 
shoved  that  the  first  difference  of  the  data  had  the  least 


amount  of  white  noise.  To  make  the  data  more  stationary,  the 
first  difference  was  accomplished  on  the  data  and  plotted  on 
FORECAST  MASTER  (Figure  2).  The  plot  now  appears  stationary 
with  a  constant  mean  and  variance. 

Now  that  the  data  are  stationary,  identification  of 
potential  models  takes  place.  Plots  and  tables  of  the  ACF 
and  PACF  for  lags  1  -  36  of  the  first  differenced  data  are 
shown  in  Figures  3  and  4  respectively.  The  ACF  and  PACF 
plot  shows  significant  spikes  at  lags  1,  3,  4,  8,  and  12. 
Significant  spikes  extend  beyond  1/2  the  standard  error  for 
lags  1-3;  one  standard  error  for  lags  4-6;  two  standard 
errors  for  lags  7  -  12;  and  three  standard  errors  for  lags 
greater  than  12  (4).  The  significant  spikes  at  4,  8,  and 
12  indicate  a  possible  seasonality  every  four  time  periods. 
The  significant  spikes  at  12,  24,  and  36  also  indicate  a 
possible  12  period  seasonality.  Lags  1  and  3  indicate 
possible  weak  moving  average  and  autoregressive  patterns. 
Based  on  the  examination  of  the  ACF  and  PACF,  possible  start¬ 
ing  models  for  the  iterative  process  of  identifying  the  most 
parsimonious  (simple)  model  were  the  ARIMA's 

Ln  (1,1,1)*(0,0,1)  , 

4 

Ln  (1,1,0)*(0,0,1)  *(0,0,1)  ,  and 

4  12 

Ln  ( 0, 1, 1 ) * ( 0, 0, 1 )  *(1,0,0)  . 

4  12 
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The  ARIMA  Ln  ( 0, 1, 1 ) * ( 0, 0, 1 )  *(1,0,0) 

4  12 

function  with  a  difference  of  order  1,  and  an  MA  o 1  order  1 
in  trc  first  set  of  parentheses.  The  second  set  of 
parentheses  indicates  an  MA  seasonal  term  of  order  4.  The 
third  set  of  parentheses  indicates  a  seasonal  AR  term  of 
order  12.  Once  possible  models  are  selected,  the  next  step 
is  to  run  the  estimations  to  see  which  model  has  the  lowest 
residual  mean  square  and  Chi-square  value. 

Flying  Hours  Estimation .  Parameter  values  for  the  ARIMA 
models : 

Ln  (1,1,1)*(0,0,1)  , 

4 

Ln  (1,1,0)*(0,0,1)  *(0,0,1)  ,  and 

4  12 

Ln  (0,1,1)*(0,0,1)  *(1,0,0) 

4  12 

were  estimated  using  the  TIMES  software.  The 

ARIMA  Ln  ( 0 , 1, 1 ) * ( 0, 1, 1 )  *(1,0,0)  had  the  lowest  residual 

4  12 

mean  square  (25.3  versus  32  and  36.7)  and  the  lowest 
chi-square  value  (25.7  with  36  degrees  of  freedom)  of  the 
three  models.  The  parameter  values  estimated  by  TIMES  for 
the  AR  and  MA  parameters  are  shown  in  Figure  5.  Once  the 
model  is  selected,  the  next  step  is  diagnostic  checking 
where  several  tests  are  performed  on  the  chosen  model  to 
test  its  validity. 
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SUMMARY  OF  MODEL  1 


HWHHHwmwmmimmwHiiHmHiHHWHHHWHwwMmHwmmwwMHmwwwHHmMWHiH 
DATA  -l  -  c135  flying  hours  1979-1985  94  OBSERVATIONS 

DIFFERENCING  ON  Z  -  1  OF  ORDER  1 

NOOEL  DEVELOPED  UITH  TRANSFORMED  DATA  :  L0G(Z(T>*  ,00000£»00) 

IIIIIUIIimillllllHIIIHIIIIHHIIlIHlIHUIIHIIIIIIIHIIimillllHmiHIHIHIIIHIIHHmlHIIIIIIIIHIH 

PARAMETER  PARAMETER  PARAMETER  ESTIMATED  95  PER  CENT 

NUMBER  TYPE  ORDER  VALUE  LOWER  LIMIT  UPPER  LIMIT 


1  AUTOREGRESSIVE  1  12  ,70994E*00  ,54041E*00  85951E*00 

2  MOVING  AVERAGE  1  1  78543E*00  .64276E.00  92909E*00 

3  MOVING  AVERAGE  2  4  .34907E*00  . 12349E*00  57445E+00 


OTHER  INFORMATION  AMO  RESULTS 

ffaHHIHMIltlfHHfMflHMflHHHflRltHHiHIIHtHIMIHIHHtlHfHHtHfmHtHMHlHHIHaHIIHHtllHI 

RESIDUAL  SUM  OF  SQUARES  20208E+00  80  D.F.  RESIDUAL  MEAN  SQUARE  .25241E-02 

NUMBER  OF  RESIDUALS  83  RESIDUAL  STANDARD  ERROR  .50240E-01 


Figure  5.  Flying  Hours  Parameter  Values/Model  Summary 
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Flying  Hours  Diagnostics ■  The  chosen  model  for 


diagnostic  checking  is  the  ARIMA  Ln  ( 0, 1, 1 ) * ( 0, 0, 1 )  * 

4 

(1,0,0)  . 

12 


The  model  equation  is: 

12  4 

(1-*B  ) ( 1-B ) ( LnX  )  =  ( 1-0B ) ( 1-0C  )e  (4) 

t  t 

where : 

B  =  the  backshift  operator  and  is  equivalent  to  X 

t-n 

and  n  is  the  superscript  above  each  B  in  the 
Box-Jenkins  notation. 

i>  =  the  seasonal  autoregressive  parameter  value 
0  =  the  moving  average  parameter  value 
0  =  the  seasonal  moving  average  parameter  value. 


Performing  the  diagnostic  checks  validates  the  model  to 
ensure  it  follows  the  pattern  of  the  data.  Some  of  the 
checks  performed  were  to  examine  the  residual  ACF  and  FACF, 
to  compare  the  chi-square  value  at  a  95%  confidence  interval, 
and  to  analyze  the  cumulative  periodogram,  histogram,  and 
power  spectrum  of  the  residuals. 

Tables  2  and  3  show  the  standard  error  values  for  the 
residual  PACF  and  ACF.  No  significant  spikes  remain,  indica- 
ing  the  model  cannot  be  significantly  improved  by  adding  more 
parameters.  The  Chi-square  test  determines  whether  the  group 
of  values  of  the  autocorrelations  is  significantly  different 
from  zero.  If  the  Chi-square  value  from  a  table  of  computed 
Chi-square  values  at  various  confidence  levels  is  larger  than 
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the  Chi-square  value  estimated  for  the  model,  then  the  data 
are  random  with  no  pattern  (3:291).  This  test  is  the  Port¬ 
manteau  Lack  of  Fit  Test.  For  the  flying  hours  ARIMA  model 
chosen,  the  Chi-square  value  is  25.7  for  36  degrees  of  free¬ 
dom.  The  table  value  for  a  95%  confidence  interval  and  36 
degrees  of  freedom  is  greater  than  43  (3:522). 

Figure  6  shows  the  Cumulative  Periodogram.  Pure  white 
noise  would  result  in  a  perfect  diagonal  line  through  the 
periodogram.  This  figure  shows  a  good  periodogram  indica¬ 
ting  just  white  noise  is  present.  Figure  7  shows  the  histo¬ 
gram  of  residuals.  It  has  a  normal  shape  and  also  confirms 
the  residuals  are  white  noise  as  desired.  Figure  8  is  the 
Power  Spectrum  which  also  supports  the  other  tests  in  indica¬ 
ting  the  residuals  are  white  noise.  A  good  power  spectrum  is 
determined  by  placing  a  ruler  horizontally  across  the  spec¬ 
trum.  If  the  ruler  picks  up  the  connecting  points  with  no 
breaks,  the  spectrum  is  acceptable  as  this  one  is. 

To  summarize,  the  ARIMA  Ln  ( 0 , 1 , 1 ) * ( 0 , 0 , 1 )  *(1,0,0) 

4  12 

is  the  univariate  model  selected  to  be  the  independent/input 
variable  used  for  identifying  the  final  transfer  function 
models.  The  next  step  is  to  prewhiten  the  dependent/output 
variable  which  is  spares  demand. 

Soares  Demand  Prewhitened .  Since  the  flying  hours 
portion  of  the  transfer  function  (with  possible  noise  para¬ 
meters)  model  is  used  to  predict  the  spares  demand  as  the 
independent  variable,  it  is  not  necessary  to  determine 
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the  actual  parameter  values  for  the  demand  data. 

Prevhitening  of  each  of  the  ten  spares  consisted  of  calcu¬ 
lating  the  first  difference  of  each  data  set  since  it  had  the 
most  stationary  data.  The  demand  data  is  combined  with 
the  flying  hours  model  and  formulation  of  the  multivariate 
transfer  models  begin  for  each  of  the  spares. 

Table  2.  Flying  Hours  PACF  of  Residuals 


PARTIAL  AUT OCORRELAT IOMS 

DATA  -  THE  ESTIMATED  RESIDUALS  -  »DEl  1 

DIFFERENCE  -  ORIGINAL  SERIES  IS  YOUR  DATA. 

DIFFERENCES  BELOW  ARE  OF  ORDER  1 

ORIGINAL  SERIES 

HEAN  OF  THE  SERIES  =-.33023E-02 
ST.  DEV.  OF  SERIES  :  .49532E-01 
NUMBER  OF  OBSERVATIONS  =  83 


1-  12 

-.02 

-.07 

-.03 

.00 

.10 

.07 

.00 

.03 

.19 

-.03 

.11 

-.15 

13-  24 

-.20 

.02 

-.04 

-.08 

.04 

-.10 

.02 

-.12 

.00 

-.17 

.04 

-.02 

25-  34 

.17 

-.10 

.03 

-.03 

.11 

-.02 

.05 

-.10 

-.14 

-.09 

.13 

.04 

37-  37  .00  .09  .01 
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Table  3.  Flying  Hours  ACF  of  Residuals 


AUTOCORRELATION  FUNCTION 

DATA  -  THE  ESTIMATED  RESIDUALS  -  NODEL  1 

OIFFERCNCINC  -  ORICINAL  SERIES  IS  VOUR  DATA. 

DIFFERENCES  BELOW  ARE  OF  ORDER  1 


ORICINAL  SERIES 

HEAH  Of  THE  SERIES  =-.33023E-02 
ST.  DEV.  OF  SERIES  =  .49532E-01 


NUNBER  OF  OISERUATIONS  = 

83 

1-  12 

-.02 

-.07 

-.02 

.01 

.10 

.07 

-.02 

.02 

.18 

-.03 

.09 

-.15 

ST.E. 

.11 

.11 

.11 

.11 

.11 

.11 

.11 

.11 

.11 

.12 

.12 

.12 

13-  24 

-.19 

.06 

.01 

-.04 

.03 

-.11 

-.01 

-.07 

-.08 

-.19 

.05 

-.03 

ST.E. 

.12 

.12 

.12 

.12 

.12 

.12 

.12 

.12 

.13 

.13 

.13 

.13 

25-  34 

.14 

-.06 

-.12 

.01 

.04 

.02 

-.03 

-.03 

-.15 

.10 

.11 

.06 

ST.E. 

.13 

.13 

.13 

.13 

.13 

.13 

.13 

.13 

.13 

.14 

.14 

.14 

37-  39 

-.09 

.04 

.04 

ST.E. 

.14 

.14 

.14 

NEAN  DIVIDED  BV  ST. 

ERROR 

r  .40740E+00 

TO  TEST  WHETHER  THIS  SERIES  IS  WHITE  NOISE,  THE  VALUE  .25731E*02 
SHOULD  IE  CONPARED  WITH  A  CHI-SQUARE  VARIABLE  WITH  3 4  DEWEES  OF  FREEDOH 
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Figure  6.  Flying  Hours  Cumulative  Periodogram 
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Figure  7.  Flying  Hours  Histogram  of  Residuals 
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POEMHITEHED  035  flyinj  hours  1978-1985 
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Figure  8.  Flying  Hours  Power  Spectrum 
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Transfer  Function  Model  Formulation .  For  an  example 
of  the  steps  taken  to  form  a  transfer  function  model  once  the 
univariate  models  are  prewhitened,  the  TAB74  spare's  data 
are  used.  The  computer  output  and  model  parameter  results 
for  the  other  nine  spares  are  located  in  Appendix  B. 

To  preview,  a  second  form  of  the  general  equation  for  a 
transfer  function  model  shown  in  Eq  (2)  is: 

r  s 

(1-s  B  .  . -s  B  ) Y  =  (w  -w  B- .  .  .-v  B  )X  +  N  (5) 
1  r  t  0  1  s  t-b  t 


where : 

r  =  the  number  of  cogent  past  Y  and  is  derived 

t 

from  the  pattern  of  the  cross  correlation.  The 
value  of  r  should  equal  the  number  of  "start-up" 
values  before  the  b  lag. 

b  =  the  number  of  cross  correlations  equal  to  zero 
before  a  significant  cross  correlation  is 
observed.  In  other  words,  b  is  the  number  of 
periods  until  the  effect  of  a  change  in 
X  (flying  hours)  affects  Y  (demand), 
t  t 

s  =  the  past  X  values.  Y  and  X  are 
t  t  t 

the  same  as  those  defined  in  Eq  (2). 


The  model  is  denoted  as  an  (r,s,b)  model  (3:377-379). 

To  identify  a  transfer  function  model,  the  steps  consist  of: 


(1)  deriving  rough  estimates  of  the  impulse  response 
weights  (Q  )  -  output  in  TIMES. 

j 

(2)  use  values  from  (1)  to  estimate  r,  s,  and  b. 

(3)  use  values  from  (1)  and  (2)  to  obtain  initial 
estimates  of  s  and  w  (3:378). 
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The  values  o£  r ,  s,  and  b  are  estimated  by  using  the  cross 
correlations  and  impulse  response  weights  from  TIMES. 

TAB74  Transfer  Function  Development .  In  developing 
the  TAB74  transfer  function  model  to  forecast,  the  flying 
hours  model  is  combined  with  the  TAB74  monthly  data  in  a 
TIMES  Fortran  program.  The  same  steps  of  identification, 
estimation  and  diagnostic  checking  are  performed. 

First,  an  identification  run  is  done  on  TIMES  to  obtain 
the  impulse  response  weights  and  cross  correlation  values. 
Table  4  shows  the  values  TIMES  calculated  for  the  impulse 
response  weights  for  TAB74.  Lag  1  has  a  large  negative  num¬ 
ber  and  this  value  appears  in  the  plot  of  the  cross  correla¬ 
tions  shown  in  Figure  9.  Table  5  shows  the  cross  correlation 

values  which  are  also  plotted  in  Figure  9.  The  series  1  data 

in  the  left  column  of  Table  5  appear  on  the  right  hand  side 

of  the  cross  correlation  plot  in  Figure  9.  It  i3  the  right 

hand  side  of  zero  which  is  used  to  help  determine  r,  s,  and 
b  along  with  the  plot  of  the  estimated  impulse  response 
weights  in  Figure  10. 


Table  4. 


Estimated  Impulse 

K 

0 

1 

2 

3 

4 

5 

6 


Response  Weights 
V1K± 

-.847 

-2.64 

.121 

-1.73 

.690 

.649 

.867 
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Figure  9.  Plot  of  the  Cross  Correlations 
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Figure  10.  Plot  of  Estimated  Impulse  Response  Weights 


Table  5 


Cross  Correlation 


Series  1  -  Prewhitened  C135  flying  hours  1978-1985 
Series  2  -  Prewhitened  demand  data  for  TAB74 


Number  of  Lags 
on  Ser  ies  1. 

0 

1 

2 

3 

4 

5 


Cross 

Correlations 

-.035 

-.109 

.005 

-.071 

.028 

.027 


Number  of  Lags 
on  Series  2 

0 

1 

2 

3 

4 

5 


Cross 

Correlations 

-  .035 
-.067 
-.010 
-.034 
-.116 
-.046 


The  large  spike  at  zero  in  Figure  9  indicates  v 

0 

is  not  zero  which  in  turn  indicates  b  =  0  for  the  TAB74 

model.  When  b  =  0,  there  is  no  time  period  delay  for  the 

effect  of  X  (flying  hours)  on  Y  .  Since  there  is  no  ap- 
t  t 

parent  pattern  in  the  right  hand  side  correlation  in  Figure  9, 

r  =  0.  For  the  model,  s  was  also  determined  to  =  0. 

Once  r,  s,  and  b  are  estimated,  the  s's  and  w's  in 

Eq  (5)  are  determined  with  a  series  of  equations  for  the 

impulse  response  functions  (3:383).  For  this  model,  only 

the  w  value  is  present  since  b  =  0.  The  value  for 
0 

w  was  estimated  to  equal  -.847  from  the  estimated 
0 

impulse  response  weights  in  Table  4.  All  other  w  values 

s 

are  zero  so  each  term  with  a  w  value  other  than  w 

s  0 

drops  out  of  Eq  (5).  The  X  terms  in  Eq  (5)  reduce  to 

t-b 
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-.867  Ln  X  .  Since  r  =  0,  the  left  hand  side  of  Eq  (5) 
t 

reduces  to  Y  . 

t 

Second,  a  TIMES  test  run  is  done  using  the  s's  and 
w's  as  input  and  output  lag  values  for  the  transfer  function 
model.  The  test  run  output  shows  whether  or  not  noise 
parameters  are  necessary.  Noise  parameters  are  needed  if 
more  than  white  noise  remains  in  the  model  (significant 
spikes  appear  in  the  noise  output).  The  TAB74  test  run 
ACF  and  PACF  had  significant  spikes  at  lags  1  and  2.  To 
prewhiten  the  noise  function  an  AR  and  MA  term  were  added  to 
the  equation  as  noise  parameters.  The  final  model  with  noise 
parameters  included,  next  goes  through  diagnostic  checking 
for  validity. 

The  diagnostic  checks  resulted  in  acceptance  of  the 
model  for  the  noise  function.  The  residual  ACF  and  PACF 
showed  no  significant  spikes.  The  Chi-square  value  of  24  at 
33  degrees  of  freedom  was  acceptable.  The  residual  mean 
square  was  low  at  .56.  Figures  11  and  12  show  both  the 
Cumulative  Periodogram  and  Histogram  of  residuals  are  good. 
The  input  lag  and  noise  parameter  values  are  shown  in 
Figure  13.  For  the  TAB74  model,  the  noise  portion  of  Eq  (5) 
becomes : 

N  =  - . 254N  +  . 157N  +  a  (6) 

t  t-1  t-2  t 

where : 

a  =  random  error  term, 
t 
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Figure  11.  Transfer  Function  Cumulative  Periodogram 
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1 

DATA  -  l  •  d33  flying  hours  1978*1985 

96  OBSERVATIONS 

Y  *  deund  data  tab74 

DIFFERENCING  ON  Y  -  1  OF  ORDER  1 

DIFFERENCING  ON  I  *  1  OF  ORDER  1 

MODEL  DEVELOPED  N1TH  TRANSFORMED  DATA  =  L0G(X(T)+  .OOOOOE+OO) 

NOISE  MODEL  PARAMETERS 

PARAfCTER 

NUMBER 

PARAMETER 

TYPE 

PARAMiTER 

ORDER 

ESTIMATED 

VALUE 

LONER  LIMIT 

95  PER  CENT 

UPPER  LIMIT 

1 

AUTOREGRESSIVE  1 

1 

-.25443E+00 

-.46367E+00 

451B1E-01 

2 

MOVING  AVERAGE  1 

2 

•  ■■■■mioatl  ■  •  »  *  i 

. 15681E+00 

-.56972E-01 

■■■■■■■■■■■■■■a 

.37059E+00 

TRANSFER  FUNCTION  PARAMETERS 

f TTITl f 1 III TTtTTTTI 

1 TT  fill TTTf f f TTTf f t«T1 

a.aMaaaaaaaiaa  a  a  aaaaaaaaa 

Mil  III  I  MTItTI  M  II 1 11 1 1 

3 

INPUT  LAG  1 

0 

-.64132E+D0 

-.24029E+01 

.11202E+01 

WHmWHWWHWHHHIHHHH  I  HHWWWtl »  WWHHHWmWHHHHWWWWWHHHWHWH 

OPTIMUM  VALUE  OF  B  IS  0 

0TOR  IffORMATION  AMD  RESULTS 


RESIDUAL  SUM  OF  SQUARES  .31011+02  91  O.F. 

NUMBER  OF  RESIDUALS  94 


RESIDUAL  JCAN  SQUARE 
RESIDUAL  STPKDMD  ERROR 


■S4057E+00 

•74871E+00 


Figure  13.  Transfer  Function  Model  Parameters 
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The  entire  equation  for  the  transfer  function  with 
parameter  values  included  for  the  TAB74  spare  is: 

Y  =  - . 867LnX  -  .615LnX  +  .61LnX 

t  t-1  t-12  t-13 

+ . 681e  +  . 303e  -  .229e 

t-1  t-4  t-5 

+  e  -  .  2  54N  +  . 157N  +  a  (7) 

t-1  t-2  t 

Each  of  the  other  nine  spares  models  were  processed 
through  the  same  steps  of  identification,  estimation,  and 
diagnostic  checking  to  develop  a  model  for  forecasting 
demand.  The  plots  and  parameter  values  for  the  remaining 
nine  spares  are  in  Appendix  B.  The  equations  for  the  nine 
models  are  listed  below: 

(1)  FAIR01 

Y  =  - . 33LnX  -  .224LnX  +  .264e 

t  t-1  t-14  t-1 

+  .132e  -  .106e  -  .32e 

t-4  t-5  t 

-  . 3N  +  . 7N  +  a 
t-x  t-6  t 

(2)  ASSY31 

Y  =  3 . 96LnX  +  2 . 69LnX  -  2.69LnX 

t  t-1  t-12  t-13 

-  3.16e  -  1.58e  +  1.27e 

t-1  t-4  t-5 

+  2 . 9  6e  +  .  1 8N  +  .  7N  +  a 

t  t-3  t-6  t 
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(3)  FRAME53 


Y  =  - . 249LnX  -  . 169LnX  +  .169LnX 

t  t-1  t-12  t-13 

+  .199e  +  .10e  -  .08e 

t-1  t-4  t-5 

-  . 249e  +  . 02N  +  a 

t  t-2  t 


(4)  TAIL84  -  NO  NOISE  PARAMETERS  NEEDED 

Y  =  -2 . OLnX  -  1 . 36LnX  +  1.36LnX 

t  t-2  t-13  t-14 

+  1.6e  +  . 8e  -  .64e  -  2e  +  a 

t-1  t-4  t-5  t  t 


(5)  COWL 8 3 

Y  =  - - 272LnX  -  .185LnX 
t  t-1  t-12 

+  . 218e  +  . le 

t-1  t-14 

-  . 127N  -  . 1 2 8N 

t-1  t-2 


(6)  SLVE85 

Y  =  . 132LnX  -  .661LnX  +  .09LnX 

t  t-1  t-2  t-12 

-  . 54LnX  +  . 4  5LnX  -  .771e 

t-13  t-14  t-1 

+  .132e  +  ,53e  -  .06e  +  .32e 

t  t-2  t-4  t-5 

-  . 21e  -  .  0  6  3N  +  . 0 57N 

t-6  t-1  t-2 

-  .  3N  +  a 
t-3  t 


-  . 185X 

t-13 

-  . 087e  -  . 27e 

t-5  t 

-  . 1 7  4N  +  a 

t-3  t 
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(7)  ASSY32 


Y  =  6 . 33LnX  +  4.3LnX  -  4 . 3LnX 

t  t-1  t-12  t-13 

-  5  .  le  -  2 . 5e  +  2e 

t-1  t-4  t-5 

+  6 . 33e  -  . 38eN  +  a 

t  n-3  t 


(8)  ASSY33  -  NO  NOISE  PARAMETERS  NEEDED 

Y  =  1 . 2LnX  +  . 82LnX  -  .82LnX 

t  t-1  t-12  t-13 

-  . 97e  -  . 48e  +  . 39e  +  1 . 2e  +  a 

t-1  t-4  t-5  t  t 


(9)  FORK 3 5 

Y  =  - . OlLnX  +  . 61LnX  +  1 . 4LnX 

t  t-4  t-6  t-7 

-  . 007LnX  +  . 4  2LnX  +  . 56LnX 

t-16  t -17  t-18 

-  . 97LnX  +  .Ole  +  .62e 

t-19  t  t-1 

+  .93e  -  l.le  -.24e  -  .37e 

t-2  t-4  t-5  t-6 

+  .  4 5e  +  . 17N  +  . 32N  +  a 

t-7  t-1  t-6  t 


These  ten  transfer  function  models  forecast  spares 
demand  for  the  Jan  -  Mar  1986  quarter.  Results  follow. 
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Mult ivar late  Forecasting  Results .  Table  6  shows  the 
Box-Jenkins  transfer  function  models  ranged  widely  in  their 
accuracy  for  predicting  quarterly  demand.  Four  of  the  models 
had  less  than  a  10%  mean  absolute  percentage  error  from  the 
actual  demand.  Each  of  the  four  (TAB74,  TAIL84,  ASSY32,  and 
ASSY33 )  overpredicted  the  demand  versus  underestimating  it. 
The  ASSY31  model  had  a  10.4%  error  and  also  overestimated. 

The  other  five  models  ranged  in  accuracy  error  from  17.6%  to 
50%  and  in  two  of  the  five  instances,  the  model  under¬ 
estimated  demand. 

The  next  models  analyzed  are  the  ten  Box-Jenkins  uni¬ 
variate  models  which  use  o.  ly  past  demand  and  the  within 
relationships  to  estimate  future  demand.  The  steps  of 
identification,  estimation,  and  diagnostic  checking  are 
shown  for  the  TAB74  model  followed  by  the  formulae  and  fore¬ 
casts  for  the  ten  spares  models. 


Table  6 


Transfer 

Function  Models' 

Forecasting 

Results 

NSN 

Actual  Value 

Forecast 

MAPE ( % ) 

FAIR01 

77 

99 

28.6 

ASSY31 

230 

254 

10 . 4 

TAB7  4 

26 

27 

3.8 

FRAME53 

34 

40 

17.6 

TAIL84 

9 

9 

0 

COWL 8 3 

10 

5 

50 

SLVE85 

13 

9 

30.8 

ASSY32 

262 

279 

6 . 5 

ASSY33 

170 

174 

2 . 4 

FORK  3 5 

17 

21 

23.5 
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Univariate  Model  Formulation 

Uni var iate  Mode  1  Identification.  Each  of  the  ten 
spares  models  were  built  using  the  same  Box-Jenkins  steps  as 
those  used  to  form  the  univariate  flying  hours  model.  The 
first  steps  in  building  the  univariate  model  TAB74  were  to 
plot  the  original  data  and  run  a  TIMES  Identification  run 
to  study  the  ACF  and  PACF  and  to  compare  the  Chi-square 
values  for  the  original  data  and  the  first  and  second  differ¬ 
ence.  The  results  indicate  whether  the  data  is  stationary 
or  not. 

The  plot  of  the  original  data  (Figure  14)  did  not  have  a 
constant  mean  indicating  the  possibility  of  non-stat ionar ity . 
The  ACF  plot  did  not  die  out  quickly  which  also  supported 
non-stat ionar * ty.  The  third  indicator  of  non-stat ionar ity 
was  the  Chi-square  value.  The  Chi-square  value  for  the  first 
differenced  data  was  the  smallest  indicating  first  differenc¬ 
ing  the  data  would  produce  the  most  stationary  data. 

The  TAB74  data  was  first  differenced.  The  ACF  and  PACF 
showed  significant  spikes  at  lags  1,  4,  6,  and  7.  The  plot 
of  the  differer.-ed  data  showed  a  much  more  stationary 
pattern  with  a  constant  mean  (Figure  15).  Based  on  the  ACF 
and  PACF,  possible  models  to  try  were  the  ARIMA's  (1,1,1), 
(0,1,2),  and  (0,1,3). 
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Unlvar late  Model  Estimation  and  Diagnostics .  A  com¬ 
parison  of  each  of  the  three  ARIMA  models  resulted  in  the 
ARIMA  (0,1,3)  having  the  smallest  Chi-square  value  of  14.0 
at  33  degrees  of  freedom  and  the  smallest  residual  mean 
square  of  0.464  (Figure  16).  To  confirm  the  ARIMA  (0,1,3) 
was  valid,  diagnostic  checks  were  done.  The  ACF  and  PACF 
each  had  no  significant  spikes  remaining.  The  Chi-square 
value  of  14  at  33  degrees  of  freedom  was  veil  within  the 
standard.  The  Cumulative  Periodogram,  Histogram  of  Residuals, 
and  Power  Spectrum  all  had  positive  results  (see  Figures 
17,  18,  and  19).  The  parameter  values  for  the  chosen 
ARIMA  (0,1,3)  are  shown  in  Figure  16.  The  equation  with 
parameter  values  is: 

Y  =  Y  +  e  -  . 349e  -  .103e  +  .425e  (8) 

t  t-1  t  t-1  t-4  t-6 

The  output  for  each  of  the  nine  Box-Jenkins  univariate 
models  is  in  Appendix  C.  The  nine  equations  with  parameter 
values  for  the  univariate  models  are  shown: 

(1)  FAIR01  -  ARIMA  (2,1,2) 

Y  =  . 804Y  -  . 106Y  -  .136Y 

t  t-1  t-2  t -  3 

+  e  -  .898e  +  . 59e 

t  t-4  t-6 
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surti«Y  of  nodel  i 


DATA  -  Z  =  demand  data  96  OBSERVATIONS 

DIFFEROCJN6  ON  2  -  1  OF  ORDER  1 

HWWWWWWWtWfWWmmWWHHHWHHWWHWlWHWHIHWHWmHWWWWHHWW 


PARAPETS? 

PARAfETER 

PARAMETER 

ESTIMATED 

95  PER  CENT 

NUHBER 

TYPE 

ORDER 

VALUE 

LONER  UNIT 

UPPER  LIMIT 

i 

ICVING  AVERAGE  1 

1 

.34861E+00 

.I6856E+00 

.52864E+CO 

** 

4 

WING  AVERAGE  1 

4 

. 1025CE+00 

-.80312E-01 

.28532E+00 

y 

N> 

WING  AVERAGE  I 

6 

-.42476E+00 

-.60660E+00 

-.24093E+00 

H<tHHW»«HWW«HHWWWWHWHHHHWWtWW>»WWtmnmWtHtW>WWWWHWHW 
OTHER  HfORHATICN  AND  RESULTS 


RESIDUAL  SUH  OF  S8UARES  .42671E+C2  92  O.F. 

NUHBER  OF  RESIDUALS  95 


RESIDUAL  l£t»  SQUARE  ,46082E*00 

RESIDUAL  STANDARD  ERROR  .68104E+00 


Figure  16.  TAB74  Parameter  Values  -  Univariate  Model 
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Figure  18.  Histogram  of  Residuals  -  TAB74  Univariate 
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Figure  19.  Power  Spectrum  -  TAB74  Univariate 
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(2)  ASSY31  -  ARIMA  (0,1,2) 

Y  =  Y  +  e  +  .197e  +  .321e 

t  t-1  t  t  -  3  t- 


(3)  FRAME53  -  ARIMA  (0,1,2) 

Y  =  Y  +•  e  -  .  146e  -  .135e 

t  t-1  t  t-1  t- 


(4)  TAIL84  -  ARIMA  (0,1,3) 

Y  =  Y  +  e  -  .077e  +  .133e 

t  t-1  t  t-1  t- 


(5)  SLVE85  -  ARIMA  (2,1,2) 

Y  =  . 56 Y  -  . 34Y  -  .329Y 

t  t-1  t-2  t-3 

+.32Y  +e  -  .  8  3e  + 
t-4  t  t-2 


(6)  COWL 8 3  -  ARIMA  (0,1,1) 

Y  =  . 89 Y  -  . 126Y  -  .116Y 

t  t-1  t-2  t-3 

-  . 1 2  6Y  +  e  +  .  4  2  7  e 

t-3  t  t-3 


(7)  ASSY32  -  ARIMA  (0,1,1) 

Y  =  Y  +  e  +  .  38  3e 
t  t-1  t  t-3 


5 


-  . 237e 
3  t-5 


.  49e 

t-6 


+  . 114e 

t-6 
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(8)  ASSY33  -  ARIMA  (0,j.,1) 

Y  =  Y  +  e  +  . 413e 

t  t-1  t  t-6 

(9)  FORK 3 5  -  ARIMA  (2,1,3) 

Y  =  Y  +  . 832Y  -  .21  Y  +  .143Y 

t  t-1  t-3  t - 4  t-5 

+  e  +  . 16  3e  -  .112e  +  .282e 

t  t-1  t-2  t-6 

Each  of  the  univariate  models  was  used  to  forecast 
demand  for  the  Jan  -  Mar  1986  quarter.  These  univariate 
models  were  developed  using  demand  relationships  to  predict 
demand.  The  previous  multivariate  models  used  the  relation¬ 
ships  among  and  between  flying  hours  and  demand  to  predict 
future  demand.  Table  7  shows  the  results  of  the  Box-Jenkins 
univariate  models'  results  for  the  ten  spares  forecasts. 

The  results  of  the  forecasts  ranged  from  three  models 
(ASSY31,  COWL83,  and  SLVE85)  predicting  the  exact  demand  while 
the  other  seven  models  ranged  in  mean  absolute  percentage 
error  from  11.1%  to  \9.b%.  Note  the  range  of  error  is 
smaller  than  the  range  for  the  transfer  function  models',  but 
then,  the  transfer  function  raoiels  had  five  with  less  than  a 
10%  error.  In  each  univariate  model,  the  forecast  was  either 
exact  or  an  overestimation.  For  aircraft  recoverable  spares, 
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Table  7 


Univariate  Models'  Forecasting  Results 
NSN  Actual  Value  Forecast  MAPE(%) 


FAIR01 

77 

92 

19.5 

ASSY31 

230 

231 

0 

TAB7  4 

26 

30 

15.4 

FRAME53 

34 

40 

17.6 

TAILS 4 

9 

10 

11.1 

COWL 8 3 

10 

10 

0 

SLVE85 

13 

13 

0 

ASSY32 

262 

293 

11.8 

ASSY33 

170 

194 

14.1 

FORK 3 5 

17 

20 

17.6 

it  is  more  mission  essential  to  overestimate  versus  under¬ 
estimate  future  requirements.  The  next  section  reviews  the 
forecasting  results  found  using  simple  exponential  smoothing 
on  FORECAST  MASTER  software. 

S imple  Exponential  Smoothing 

The  *-en  spares  demand  data  sets  were  run  on  FORECAST 
MASTER  software  using  an  a  =  0.4  value  for  each  run.  This 
value  was  recommended  by  Craig  Sherbrooke  in  his  study 
discussed  in  Chapter  II.  The  results  of  the  forecasts  for 
the  Jan  -  Mar  1986  quarter  are  in  Table  8.  Note  that  three 
models  (TAIL84,  ASSY3-,  and  ASSY33)  forecast  exact  demand. 
Five  of  the  other  models  were  within  10%  of  the  actual  value. 
The  other  two  models  (FRAME53  and  FORK35)  had  errors  cf  14.7% 
and  41.2%;  both  were  overest imat ions  of  demand. 
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Table  8 


Simple  Exponential  Smoothing  Models'  Forecasting  Results 
NSN  Actual  Value  Forecast  MAPE ( % ) 


FAIR01 

77 

81 

5.2 

ASSY31 

230 

212 

7.8 

TAB74 

26 

27 

3.8 

FRAME53 

34 

39 

14 . 7 

TAIL84 

9 

9 

0 

COWL 8 3 

10 

11 

10 

SLVE85 

13 

12 

7.7 

ASSY32 

262 

263 

0 

ASSY33 

170 

170 

0 

FORK 3 5 

17 

24 

41.2 

Answers  for  Investigative  Questions 

In  Chapter  I,  questions  were  asked  to  aid  in  the  analysis 
and  development  of  time  series  models  to  forecast  C-135 
recoverable  spares.  Below  are  the  answers  and  findings  for 
the  questions. 

(1)  What  type  of  relationship  exists  between  C-135  flying 
hours  and  demand  in  the  transfer  function  models? 

In  each  of  the  ten  multivariate  transfer  function 
models,  the  cross  correlation  values  were  not  as  high  as  ex¬ 
pected.  Low  cross  correlation  values  indicate  the  relation¬ 
ship  between  flying  hours  and  demand  for  the  C-135  data  is 
not  strong.  Two  conditions  which  are  necessary  for  a 
transfer  function  to  work  are,  the  impulse  response  weights 
must  converge  and  second,  some  small  incremantal  change  in 

X  (flying  hours)  results  in  a  small  incremental  change 
t 

in  Y  (demand )  ( 4  )  . 
t 
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In  the  majority  o£  the  transfer  function  models,  the 
spares  demand  was  related  to  flying  hours  1,  12,  and  13  time 
periods  prior.  The  transfer  function  equations  shown  earlier 
in  thi3  chapter  show  the  actual  prior  time  periods  in  which 
demand  is  related  to  flying  hours  for  each  individual  model. 
When  the  flying  hours  data  was  prewhitened,  there  was  no  trend 
in  the  data  but  a  seasonal  component  existed  at  time  periods 
of  4  and  12. 

(2)  and  (3)  Which  forecasting  techniques  best  modelled  the 
data  and  produced  the  most  accurate  results? 

Each  of  the  three  types  of  models  (Box-Jenkins  multi¬ 
variate,  Box-Jenkins  univariate,  and  simple  exponential 
smoothing),  modelled  the  data  accurately  since  each  of  the 
three  types  of  models  forecast  one  quarter  of  demand  for  the 
spares  with  accurate  results.  The  objective  in  developing  the 
Box-Jenkins  models  was  to  build  models  which  "fit"  the  actual 
data  versus  building  to  obtain  "forecast"  accuracy.  The  final 
section  compares  the  mean  absolute  percentage  error  for  each 
of  the  ten  spares  using  the  three  types  of  models. 

Table  9  is  a  comparison  of  the  mean  absolute  percentage 
errors  for  the  ten  spares.  The  Box-Jenkins  transfer  ft’-'  a 
and  univariate  demand  models  did  well  but  the  exponential 
smoothing  with  a  =  0.4  had  the  lowest  average  percentage 
error  for  the  ten  different  spares. 
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Table  9 


Comparison  of  Forecasting  Results 


NSN 

Transfer  Func . 
MAPE ( % ) 

Univar iate 
MAPE ( % ) 

Simple  Exp 
MAPE ( % ) 

FAIR01 

28 . 6 

19 . 5 

5.2 

ASSY31 

10.4 

0 

7.8 

TAB  7  4 

3.8 

15.4 

3.8 

FRAME 5 3 

17.6 

17.6 

14.7 

TAIL84 

0 

11.1 

0 

COWL83 

50 

0 

10 

SLVE85 

30.8 

0 

7.7 

ASSY32 

6 . 5 

11.8 

0 

ASSY33 

2.4 

14.1 

0 

FORK 3 5 

23 . 5 

17.6 

41 . 2 

AVERAGE  MAPE: 

17.4 

11.0 

9.0 

Conclusions  and  Recommendations 


For  the  last  several  years,  emphasis  has  been  placed  on 
finding  forecasting  techniques  to  improve  on  the  past  under¬ 
estimations  of  aircraft  recoverable  spares.  Two  views  of 
approaching  forecasting  spares  which  were  discussed  in  the 
Chapter  II  literature  review  are  to  use  time  series  techniques 
or,  techniques  such  as  exponential  smoothing.  The  purpose  of 
this  research  was  to  compare  Box-Jenkins  time  series  tech¬ 
niques  to  simple  exponential  smoothing  with  C-135  demand  and 
flying  hours  data.  The  Box-Jenkins  multivariate/transfer 
function  models  use  the  relationship  of  flying  hours  and 
demand  to  forecast  future  demand.  The  Box-Jenkins  uni¬ 
variate  models  used  past  demand  and  it’s  relationship  with 
time  and  past  demand  to  forecast  future  demand.  Simple 
exponential  smoothing  used  past  demand  data  with  more 
emphasis  placed  on  recent  demand  to  forecast  future  demand. 
Chapter  IV  provided  a  comparison  of  the  results  of  the  fore¬ 
casts  for  each  of  the  three  methods. 

Findings 

Key  findings  of  the  research  are  summarized  below: 

Low  cross  correlation  values  during  transfer  function 
model  development  indicated  flying  hours  were  not  as 
highly  correlated  with  demand  as  expected. 

The  average  mean  absolute  percentage  error  for  the  ten 
spares  was  17.4%  when  using  the  transfer  function  models. 

Two  of  the  ten  spares'  demand  were  underestimated  when 
using  the  transfer  function  models. 
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The  average  mean  absolute  percentage  error  for  the  ten 
spares  was  11%  for  the  Box-Jenkins  univariate  models. 

The  univariate  models  did  not  underestimate  demand. 

The  average  mean  absolute  percentage  error  for  the  ten 
spares  was  9%  using  simple  exponential  smoothing. 

The  simple  exponential  smoothing  models  underestimated 
demand  for  two  spares'  forecasts. 

Research  Limitations 

Box-Jenkins  time  series  analysis  requires  at  least  forty 
and  preferably  100  data  values.  Because  of  this,  only  one 
quarter's  demand  was  withheld  to  take  advantage  of  Box- 
Jenkins*  powerful  techniques  of  developing  models  which 
predict  accurate  mid-to-long  term  forecasts.  Simple  expo¬ 
nential  smoothing  is  advantageous  for  short  term  (1-3  months) 
forecasting.  Therefore,  only  having  one  quarter  of  data  to 
compare  actual  and  forecast  values  favored  exponential 
smoothing . 

Conclusions 

This  research  determined  Box-Jenkins  techniques  can  be 
used  to  accurately  forecast  future  aircraft  recoverable 
spares.  Previous  theses  built  transfer  function  models  but 
did  not  use  the  models  to  forecast.  Numerical  results  show 
the  transfer  function  models  which  use  flying  hours  to  fore¬ 
cast  demand,  may  be  more  complex  than  needed  since  both  the 
Box-Jenkins  univariate  and  simple  exponential  smoothing 
models  produced  better  results  (11%  and  9%  versus  17%)  for 
one  quarter  when  only  using  past  demand.  One  reason  for  the 
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transfer  function  models  not  forecasting  as  well  is  the 
C-135  flying  hours  were  not  as  highly  correlated  to  demand 
as  expected.  It  should  be  noted  though,  that  Box-Jenkins 
techniques  are  used  for  mid-to-long  term  (1-2  years)  fore¬ 
casting  while  simple  exponential  smoothing  is  short  term 
(1-3  months).  Due  to  data  constraints,  only  one  quarter  of 
demand  was  forecast  and  compared  to  actual;  this  favors 
simple  exponential  smoothing.  The  accuracy  of  simple  expo¬ 
nential  smoothing  with  a  =  0.4,  as  Sherbrooke  recommended, 
shows  it  may  not  be  necessary  to  remove  "within"  and 
"between"  relationships  in  the  data  to  forecast  accurately. 

The  results  of  these  forecasting  models  impact  the 
AF/DOD  budget  process  since  the  determination  of  spares  re¬ 
quirements  is  usually  made  2-3  years  ahead  of  actual  sup¬ 
port  arriving.  Box-Jenkins  techniques  are  powerful  enough 
to  forecast  1-2  years  ahead  for  demand.  In  cases  where 
item  managers  need  short  terra  forecasts,  this  research 
indicates  simple  exponential  smoothing  would  be  very 
effective.  Depending  on  the  time  horizon  of  the  forecast 
needed  and  the  amount  of  data  available,  both  the  TIMES 
and  FORECAST  MASTER  software  are  available  for  Air  Force 
use.  TIMES  was  originally  designed  for  A’r  Force  use  on 
mainframe  computers  which  compile  Fortran.  It  has  recently 
been  converted  for  personal  computer  use  at  AFIT. 

FORECAST  MASTER  software  is  personal  computer  software 
which  is  very  user  friendly  and  capable  of  doing  both 
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Box-Jenklns  and  simple  exponential  smoothing  models  as  well 
as  many  others. 

One  other  conclusion  of  this  research  is  that  Air  Force 
data  collection  and  storage  is  not  condusive  to  using 
Box-Jenkins  analysis  since  a  large  number  of  data  points  are 
needed  and  the  majority  of  data  is  stored  as  quarterly  which 
gives  four  points  per  year  of  data.  Capabilities  of  today's 
personal  computers  and  mainframes  allow  storage  of  vast 
amounts  of  data  which  fifteen  years  ago  wouldn't  fit  on  a 
mainframe  computer.  Whenever  possible,  data  should  be  stored 
at  a  minimum  of  monthly,  which  opens  the  doors  to  powerful 
techniques  such  as  Box-Jenkins.  Incomplete  data  storage  also 
lessens  the  capability  of  getting  true  and  accurate  results 
based  on  past  data  if  there  are  excessive  gaps  in  the  data. 

Recommendations 

The  models  used  in  this  thesis  produced  accurate  results, 
so  additional  research  should  be  accomplished  to  compare  the 
mld-to-long  term  forecasting  of  the  three  techniques  since 
only  short  term  was  tested.  In  studying  Box-Jenklns  time 
series  techniques,  it  is  also  recommended  that: 

-  The  f lying-hours-demand  relationship  be  studied  for 
other  weapon  systems  and  their  spares. 

-  Multivariate  models  be  further  tested  which  also 
take  into  account  demand  relationships  with  variables 
such  as  number  of  take-offs  and/or  landings,  number 
of  sorties,  and  age  of  the  aircraft  fleet. 

-  Studies  be  accomplished  which  explore  combining  spares 
based  on  similar  stock  class  or  end  item  use  to  test 
whether  or  not  aggregate  forecasting  is  possible. 
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Appendix  A:  Flying  Hour  and  Spares  Demand  Data 

Table  A1 

Monthly  Flying  Hours  Data 


Date 

Flying 

Hours 

Date 

Flying 

Hours 

Date 

7801 

21045 

8108 

20516 

8503 

7802 

21718 

8109 

23877 

8504 

7803 

24259 

8110 

23592 

8505 

7804 

22467 

8111 

21251 

8506 

7805 

23068 

8112 

19078 

8507 

7806 

22734 

8201 

19643 

8508 

7807 

21738 

8202 

20939 

8509 

7808 

24310 

8203 

23821 

8510 

7809 

24361 

8204 

22718 

8511 

7810 

24623 

8205 

21688 

8512 

7811 

22772 

8206 

22932 

7812 

18724 

8207 

22430 

7901 

22055 

8208 

21308 

7902 

21126 

8209 

20719 

7903 

24432 

8210 

22993 

7904 

23172 

8211 

21484 

7905 

23995 

8212 

19332 

7906 

22588 

8301 

20645 

7907 

21846 

8302 

20358 

7908 

24984 

8303 

23593 

7909 

21511 

8304 

22420 

7910 

24319 

8305 

21834 

7911 

21817 

8306 

22266 

7912 

17587 

8307 

20418 

8001 

21066 

8308 

22498 

8002 

20539 

8309 

20978 

8003 

21454 

8310 

24294 

8004 

22827 

8311 

20942 

8005 

22466 

8312 

18531 

8006 

21855 

8401 

21017 

8007 

22179 

8402 

20516 

8008 

21009 

8403 

22794 

8009 

19514 

8404 

23141 

8010 

23817 

8405 

22251 

8011 

20728 

8406 

20903 

8012 

19307 

8407 

21345 

8101 

20981 

8408 

23393 

8102 

21189 

8409 

21605 

8103 

22703 

8410 

24985 

8104 

22645 

8411 

21308 

8105 

21340 

8412 

17854 

8106 

21132 

8501 

21379 

8107 

21298 

8502 

19842 

Flying 

Hours 

22408 

22578 

22338 

21451 

21991 

22759 

21499 

25407 

20361 

18895 


79 


Table  A2 


Quarterly  Demand  Data 


Date 

FAIR01 

ASSY31 

TAB74  FRAME53 

TAIL84 

7803 

23 

40 

11 

12 

21 

7806 

17 

40 

13 

12 

22 

7809 

19 

43 

10 

10 

24 

7812 

17 

29 

14 

10 

44 

*7903 

16 

33 

13 

10 

42 

*7906 

17 

38 

14 

11 

29 

7909 

18 

42 

13 

11 

23 

*7912 

17 

41 

14 

10 

30 

8003 

14 

41 

15 

9 

46 

8006 

18 

46 

16 

13 

43 

*8009 

19 

49 

15 

13 

27 

8012 

21 

52 

15 

12 

26 

8103 

17 

78 

11 

20 

9 

*8106 

19 

85 

11 

21 

9 

8109 

21 

99 

10 

22 

9 

8112 

20 

100 

10 

25 

7 

8203 

18 

107 

10 

31 

6 

8206 

17 

122 

12 

32 

7 

8209 

14 

124 

9 

32 

6 

8212 

21 

138 

9 

34 

7 

8303 

28 

146 

9 

35 

4 

8306 

28 

137 

8 

45 

3 

8309 

23 

112 

10 

45 

3 

*8312 

22 

110 

10 

44 

4 

*8403 

22 

110 

10 

43 

4 

8406 

23 

109 

10 

42 

5 

8409 

28 

107 

10 

43 

6 

8412 

28 

97 

11 

38 

7 

8503 

58 

145 

20 

43 

8 

8506 

52 

154 

15 

50 

9 

8509 

85 

196 

26 

39 

9 

8512 

81 

219 

27 

37 

9 

**8603 

77 

230 

26 

34 

9 

ignates  quarters 
Ignated  quarter 

with  mi 
vithhe Id 

ssing  data 
to  compare 

with 

forecas 
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Table  A2  (CON'T) 


Quarterly  Demand  Data 


Date 

COWL  8  3 

SLVE85 

ASSY32 

ASSY33 

PORK  3 5 

7803 

53 

9 

39 

44 

42 

7806 

48 

9 

43 

44 

44 

7809 

40 

10 

52 

52 

43 

7812 

48 

30 

36 

40 

44 

*7903 

43 

27 

40 

44 

44 

*7906 

38 

25 

45 

46 

47 

7909 

34 

22 

52 

51 

47 

*7912 

40 

33 

50 

50 

44 

8003 

44 

44 

53 

52 

41 

8006 

43 

46 

64 

59 

43 

*8009 

40 

40 

70 

62 

42 

8012 

32 

35 

74 

66 

41 

8103 

27 

24 

94 

83 

45 

*8106 

21 

18 

109 

90 

43 

8109 

18 

13 

121 

105 

42 

8112 

17 

12 

117 

111 

40 

8203 

15 

13 

120 

117 

40 

8206 

15 

11 

134 

142 

43 

8209 

13 

13 

133 

144 

42 

8212 

11 

13 

144 

155 

40 

8303 

13 

13 

154 

166 

40 

8306 

15 

16 

141 

158 

33 

8309 

13 

16 

113 

131 

*7 

*8312 

15 

15 

3  12 

128 

45 

*8403 

15 

15 

711 

128 

45 

8406 

16 

14 

109 

126 

48 

8409 

16 

12 

94 

98 

43 

8412 

15 

12 

88 

88 

42 

8503 

13 

14 

124 

11? 

37 

a  unc 

W  V  w 

12 

U 

130 

40 

8509 

12 

11 

219 

178 

25 

8512 

11 

12 

280 

171 

22 

**8603 

10 

13 

262 

170 

17 

*  designates  quarters  with  missing  data 

**  designated  quarter  withheld  to  compare  with  forecasts 
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Appendix  B:  Transfer  Function  Models 

DATA  ■  I  =  cl35  flying  houri  1978-1885  86  0BSER7AT10IS 

T  1  deaand  data 

D1FFE££*CIIG  01  T  -  1  OF  ORDER  1 
DIFFEREMCIMG  OK  I  -  1  OF  ORDER  1 

MODEL  DEVELOPED  KITH  T RAISFORKED  DATA  =  LOG (X(T) ♦  OOOOOE+OO) 

«in»iititiiuit«iut«tttntitini»tttttMtinniiitttiinttnitittttnntittnnimniimutnnttiiiiittiintti 
KITE  MODEL  PARAMETERS 

lifil»llltlttilltimtlil»ltlimitttltm»tllti»lt»ltllliltltllllltllillllHlf  Itttttttlllitttttltttltllltliilltltt 

PARAMETER  PARAlffTER  PARAMETER  ESTIMATED  95  PER  CEIT 
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Figure  B1 .  FAIROl  Transfer  Function  Parameter  Values 
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Figure  B3. 


FAIROl  Transfer  Function  Histogram  of  Residuals 
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Figure  B4.  ASSY31  Transfer  Function  Parameter  Values 
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Figure  B6.  ASSY31  Transfer  Function  Histogram  of  Residuals 
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Figure  B7 .  FRAME53  Transfer  Function  Parameter  Values 
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Figure  B9 .  FRAME53  Transfer  Function  Histogram  of  Residuals 
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Figure  BIO.  TAIL84  Transfer  Function  Parameter  Values 
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Figure  Bll.  TAIL84  Transfer  Function  Cumulative  Periodograra 
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Figure  B12.  TAIL84  Transfer  Function  Histogram  of  Residuals 
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Figure  B13.  COWL83  Transfer  Function  Parameter  Values 
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Figure  B14.  COWL83  Transfer  Function  Cumulative  Periodogram 
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Figure  B15.  COWL83  Transfer  Function  Histogram  of  Residuals 
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Figure  B16.  SLVE85  Transfer  Function  Parameter  Values 
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Figure  B17.  SLVE85  Transfer  Function  Cumulative  Perlodogram 
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Figure  B18.  SLVE85  Transfer  Function  Histogram  of  Residuals 
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NOISE  MODEL  PARAMETERS 

HHlWHHIHHIHHHWHHWWHWWMHHHHWHWHHHHHWWHHi  .uni’  llHIHIIHHWIUW 

PARAMETER  PARAMETER  PARAMETER  ESTIMATED  95  PER  CENT 

NUMBER  TYPE  ORDER  VALUE  LONER  LIMIT  UPPER  LIMIT 


MOVING  AVERAGE  1  3  -.J7448E+00  -.57047E+00  -.18249E+00 


TRANSFER  FUNCTION  PARAMETERS 

HHHUWWtmHimHIUW»WtWWWUHHH)HmHHWI)IHmHmHHUmiU)HUmilimHHW« 

2  Item  LAB  l  0  .63294E+01  -.15040E+01  . 14163E+02 

OPTIMUM  VALUE  OF  8  IS  0 

WHmHHHHHUHmwUWmuWtmwtWWHWHHHHHHHHHWHHWmHIIlmWUHIHmW 
OTMER  INFORMATION  AND  RESULTS 


RESIDUAL  SUM  OF  SOUARES  -14134E+04  93  D.F.  RESIDUAL  MEAN  SQUARE  .152006*02 

NUMBER  OF  RESIDUALS  95  RESIDUAL  STANDARD  ERROR  .38988E+01 


Figure  B19.  ASSY 32  Transfer  Function  Parameter  Values 
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Figure  B21.  ASSY32  Transfer  Function  Histogram  of  Residua 
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DATA  -  X  -  cl35  Hying  hours  1978-1985 
v  =  demand  data  assy33 


96  OBSERVATIONS 


DIFFERENCING  ON  Y  -  1  OF  ORDER  1 
D1FFERENC1N6  ON  X  -  1  OF  ORDER  1 

MODEL  DEVELOPED  KITH  TRMSFORJCD  DATA  =  LOGIX(T)*  .OOOOOE+OO) 


NOISE  MODEL  PARAfETERS 

IWHItHHWIHHIWtHWHWUWHHHmtllHWHtiUWmHHHWWWIlWHWWHWHttHWWt 

PARAfETER  PARAMETER  PfWA/ETER  ESTIMATED  95  PER  CENT 

NUMBER  TYPE  ORDER  VALUE  LONER  LIMIT  UPPER  LIMIT 


wwwHiHWHHHwwHWHuwHUwimwmmiwuwmHni;:;' !  ’iih» :  i»h  ;  i  !«hkh«hwh!  ;m 
TRANSFER  FUNCTION  PARAftTERS 

HHIHHHHHHHHWHHWWHHmiHWHtHHWWHHWHHHHIWHUHWHWWHtWWWtHHH 


IfFUT  LAG  1 


.  12178E+01 


-■57504E+01 


.81861E+01 


OPTIMUM  VALUE  OF  B  IS  0 


0T1CR  IftORMATION  M)  RESULTS 

lHIIHHIUIHIItWHWtllHHmWWWWHWHWHHI«mHHHHWWWWHWH«mmHIIHUHIl 

RESIDUAL  SUM  OF  SQUARES  .98672E+03  9*  D.F.  RESIDUAL  IEAN  S8UARE  .10497E+D2 

NUMBER  OF  RESIDUALS  95  RESIDUAL  STANDARD  ERROR  .32399E+01 


Figure  B22.  ASSY33  Transfer  Function  Parameter  Values 
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Figure  B23.  ASSY33  Transfer  Function  Cumulative  Periodogram 
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Figure  B24.  ASSY33  Transfer  Function  Histogram  of  Residuals 
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DATA  -  I  =  c 135  Flyinq  hours  1978-1985 
Y  =  deund  diU  Tork35 


96  OBSERVATIONS 


DIFFERENCING  ON  Y  -  1  OF  ORDER  1 
DIFFERENCIN6  ON  X  -  1  OF  ORDER  1 

MODEL  DEVELOPED  fc'TH  TRANSFORM  DATA  =  LOG(XtT)*  .OOOOOE+OO) 


NOISE  NOOEL  PARAfCTERS 


PARAMETER 

NUMBER 


PARAMETER 

TYPE 


2 


PtfWCTER 

ORDER 


ESTIMATED 

VALUE 


95  PER  CENT 


LONER  LIMIT 


AUTOREGRESSIVE  1 
MOVING  AVERAGE  1 


.32033E+00 
.  16500E+00 


. A8745E— 01 
-.59507E-01 


UPPER  LIMIT 


«»iniiii«nniim»«minm»nHimmiimilll|HHIll|i|l|HmH»HlHHm«li  1 1  n  gg  II  I  n  1 1  ■  ■  1  M  1  ■ 

•  W  ITtlfllll  Ww  a  I  IVTI  llllllllllllllllllllirillliii  ITTTTTTITI  Illllllll  III  I  IITT1TTI  MTTTTf  1 1  ■  TIT  I  ITTwTTTTTf  r  ■  ■  IfTTTTTTf 


.59192E+00 

-3895CE+00 


HIHWHWWHIHHHWHWtHWHHWHWWWHHWWHHWHIWHUWHWWHWSSWtWmwHHH 
TRANSFER  FUNCTION  PARAfCTERS 


I  (CUT  LAG  1 
I  (CUT  LAG  1 
[(PUT  LAG  1 


12686E+01 
.61395E+00 
.  14256E+01 


-.40913E+01 

-.24714E+01 

13788E+01 


.  15541E+01 
.36993E+01 
.42301E+01 


OPTIMUM  VALUE  OF  B  IS  ♦ 

OTHER  I  (FORMATION  AND  RESULTS 

IHtHHWHWHHmiWWWWtiHMWtHttWmmWHHWIWWHHtmWHWHHtWWHHHWm 


RESIDUAL  SUM  OF  S8UARES  .92377E+02  78  D.F.  RESIDUAL  (CAN  SOUARE  .11843+01 

NUMBER  OF  RESIDUALS  83  RESIDUAL  STANDARD  ERROR  .10883+01 

Figure  B25.  FORK35  Transfer  Function  Parameter  Values 
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Figure  B26.  FORK35  Transfer  Function  Cumulative  Periodogram 
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Figure  B27.  FORK35  Transfer  Function  Histogram  of  Residuals 
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Appendix  C:  Univariate  Models 


smaUBT  OF  MODEL  1 

DATA  -  Z  •  demand  data  88  0BSES7ATI0IS 

DIFFEREMCIMG  01  Z  -  1  OF  0&DES  1 


it  iminiif  linn  in  u  iiiiiiiii  iiiiiiuim  in  iiiiimiiiiiiittiiiiiiiiiiiti  it  iiititiitiiitii  iiiiiiiiiii  iiiHitii 


parameter 

PARAMETER 

PARAMETER 

ESTIMATED 

95  PER  CEMT 

HUGER 

TTPE 

ORDER 

VALUE 

LOMER  LIMIT 

UPPER  LIMIT 

iiiiiiiiiiiiiiimiiiiiiiiiiiiHiKiiiimiiiiiiiiiiMiimiiiiiiimiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiimiiimiM 

1 

AUTOREGRESSIVE  1 

1 

-  10630EH50 

-.31513E*00 

10252E*00 

2 

AUTOREGRESSIVE  I 

3 

•  - I3600E*00 

-  3460SE'00 

74093E-0I 

3 

MOVING  AVERAGE  1 

4 

89832E-0I 

-  88349E-01 

2690 1E»00 

4 

OVIK  AVERAGE  1 

6 

-  58958E*00 

-  78853EMJ0 

•  4 !  064E*00 

IIHIIIIIIMIIIIIIIIIIIIIIIIIIIIMIIIIIIIIIIIIIIMIIIIIIUIIIIIIIIIIIIIIIIIKIIIIIIIUIIIIIIIIIIIUIIIIIIIUIIimi 

OTHER  IIF0RMATI0I  AID  RESULTS 

iiftiuiitif  itiiiiiiiiiii  ii  ■■■1111111111  if  itittiiniii  util  i  ii  mi  hi  iiuiiif  i  in  iiiiim  mu  111111111111111111  in 


RESIDUAL  SOU  OF  SQUARES  ,20158E*03  88  OF 


RESIDUAL  MEAI  SQUARE 


22906E  *0 1 


IURBER  OF  RESIDUALS  92 


RESIDUAL  STARDARD  EEROS  15I35E-01 


Figure  Cl.  FAIROl  Univariate  Parameter  Values 
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Figure  C3.  FAIROl  Univariate  Histogram  of  Residuals 
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PRENHITENED  dwand  data 
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Figure  C4.  FAIROl  Univariate  Power  Spectrum 
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SUMIART  OF  MODEL  1 

lilllllliilllHHIIIIIIIIIIIIIIIIIIfIMIIIIIIIIIIIIHillllllllllllllllllHIIIIIIIIimillllllllllHIIIIimilllllll 

DATA  -  Z  »  deiand  data  uiy31  9«  OBSEB7ATIOIS 

DIFFEBEICIKJ  01  Z  -  1  OF  OBOES  1 

lllltll  |  tilt  Hill  t||ttllt|||i«!ltttttlltltttlttltltltllllf*llll!IIIIIMItHlftltti!tl!t!i!l!tHII  lllllllllttf  III  III 

PABAMETEB  PABAMETEB  PABAMETEB  ESTIMATED  95  PEI  CE1T 

IUMBEB  TTPE  OBOES  VALUE  LOMEB  LIMIT  OPPEB  LIMIT 

nniniintitttnntntinnntnntitttnttniiutninntnntttutninntiinntnnnntuittntntitiiiintt 

1  MOVIK  A7EBA3E  1  3  -,196B2E*00  -,39780E*00  4I662E-02 

2  MDFIK  A7EBAGE  1  «  -,32083E»00  -,52257E*00  -  11909E*00 

nnninnnnininniininininininnnnnnnnnnniinnnininnniniininnniiiiiinniininnu 
OTHEB  IIFOBMATIOI  AID  RESULTS 


nninnnnnnnnnnnnnnnniiiniiinnninininninniiiiiiiniiiiiiiniiiiiiiniiiiiiiiiiiiininni 


RESIDUAL  SOM  OF  SQUARES 

. 81632E+03 

93  D.F. 

RESIDUAL  MEAI  SQOABE 

60271E+01 

IUMBEB  OF  RESIDUALS 

95 

RESIDUAL  STAID ABD  EBBOB 

,25743E*01 

Figure  C5.  ASSY31  Univariate  Parameter  Values 
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Figure  C6.  ASSY 31  Univariate  Cumulative  Periodogram 
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Figure  C7 .  ASSY31  Univariate  Histogram  of  Residuals 
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Figure  C8 .  ASSY31  Univariate  Power  Spectri 
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SUMMARY  OF  MODEL  1 

HHOtlHHWHWWHHftlltilHtHHtHHtlWlHmtHHtUHHHHWitmHWHtKHHHItiWiliMWfH 


DATA  -l  •  dewnd  data  96  OBSERVATIONS 

DIFFERENCING  ON  l  -  1  OF  ORDER  1 

WHHHHWWlHIHHWimittHHHKHtttiHHimiiillltHHtttHIHtHfHtHHIKKKHWttimwfWfH 


PARAMETER 

PARAMETER 

PARAMETER 

ESTIMATED 

95  PER  CENT 

NUMBER 

TYPE 

ORDER 

VALUE 

LONER  LIMIT 

UPPER  LIMIT 

4HHH4#HHWHHHHWWHHWm4HHWWHHWH«HWHHWHHWWHWH4WWHHWHHmmm  | 

1 

I 

MOVING  AVERAGE  1 

1 

. 14553E+00 

57480E— 01 

.34854E+00 

2 

MOVING  AVERAGE  1 

C 

J 

.13526E+00 

-.69389E-01 

.33 99 2E +00 

OIHER  INFOKJImTION  AND  RESULTS 

IHttiHitHIKHWWHWtWHWmiimttfHIIWiHttlHMiiiHtHtHtWimilHtWIWHHtWMMItHW 


RESIDUAL  SUM  OF  SQUARES  .73859E+02  93  D.F.  RESIDUE.  MEAN  SQUARE  .79413E+00 

NUMBER  OF  RESIDUALS  95  RESIDUAL  STANDARD  ERROR  .89117E+OC 


Figure  C9 .  FRAME53  Univariate  Parameter  Values 
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Figure  CIO.  FRAME53  Univariate  Cumulative  Perlodogram 
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Figure  Cll.  FRAME53  Univariate  Histogram  of  Residuals 
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Figure  C12.  FRAME53  Univariate  Power  Spectrum 
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suimy  of  model  i 

DATA  -  l  -  demand  data  9fc  OBSERVATION'S 

DIFFERENCE  ON  Z  -  1  OF  ORDER  1 

PARAttTER  PARATETER  PARAMETER  ESTIMATED  95  PER  CENT 

NUMBER  TYPE  ORDER  VALUE  USHER  LIMIT  UPPER  LIMIT 

HWHWWHHHHHtHHWHWWHWtHtmWWHtWHWtfWWHHWHtWIHHWtmtfWIWHm 

1  MOVING  AVERAGE  1  1  .77157E-CI  -.1253«E+00  .2'%cJi*00 

2  MOVING  AVERAGE  1  3  -.iC^T+OO  -.33693+00  .71B72E-01 

3  MOVING  AVERAGE  I  6  .23670E+00  .32123E-01  .44127E+00 

HHmiWIHHHIIIKIHItHHmHHilltHWtHmtWttWWtltWHIttHHWHHHtHItmtntHHtttttXH 

OTHER  INFORMATION  AND  RESULTS 

RESIDUAL  SUM  OF  SOUARES  .17554E+03  92  D.F.  RESIDUAL  *AN  SQUARE  . l<»O8Sc+0J 

NUMBER  OF  RESIDUALS  95  RESIDUAL  STANDARD  ERROR  .13813E+01 


Figure  C13.  TAIL84  Univariate  Parameter  Values 
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Figure  C14.  TAIL84  Univariate  Cumulative  Periodogram 
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Figure  C15.  TAIL84  Univariate  Histogram  of  Residuals 
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Figure  C16 .  TAIL84  Univariate  Power  Spectrum 
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SUMMARY  OF  MODEL  1 


Ifttlfllltltilltlttllltlllllttttlllllltlllttlllllltltltltlf Itlllllltllllltltllllliltiltlltltltttlttlllttllltllltltl 

DATA  -  Z  -  deaand  data  cowl  83  98  OBSEBVATIOIS 

DIFFEBEICIIG  01  Z  -  1  OF  ORDER  1 

fflf Iff  If ff Ifltlfflfff fltlltftf Itlllltlf Ittltltltlllftlltff Mllltffflf I Itllltttf fltlHflf fffllflllfllfllftltllllllt 


parameter 

PARAMETER 

PARAMETER 

ESTIMATED 

95  PER  CEIT 

IUMBEB 

TYPE 

ORDER 

VALUE 

LONER  LIMIT 

UPPER  LIMIT 

ll»»*lt»ttlllll»t«ltf||*llttlttf lilt llltllllltlftltlll lilt ttllttltl«t«*l 

[tl»lllltll»*ll*«fttlVIIIIIIIIII«ll9llflll«l 

l 

AUTOREGRESSIVE  1 

1 

*. 11558E*00 

- . 32783E+00 

96712E-01 

2 

AUTOREGRESSIVE  1 

2 

- . 12574E*00 

-  33782E+00 

8643BE-01 

3 

MOVIK  AVERAGE  1 

3 

-.212!8EM>0 

-,42B84EfOO 

24768E-02 

4 

MD7I1G  AVERAGE  I 

8 

- . 1 1390E+00 

-,32563Ef00 

97823E-01 

tfllttff  If  flflltlff  tllllf  llflltltlllf  ■Itlllllflf  llfflltltlf  fllflf  tlfftltlflf  If  llltlltlllllf  Ilf  llllllllllllllllltlll 

OTHER  IIF0BMATI0I  AID  RESULTS 

ifiiiiiiiiiiiiiiiiiiiiiiffiifiifiiiiiiiiiiiiifiiiiiiiiiiiiiiiiiiifiiiiiiiiiiifiifiiiiiiitiiiiiiiiiiifiiififiiiiiiii 


RESIDUAL  SUM  OF  SQUARES 

,70903E*02 

89  D.F. 

RESIDUAL  MEAI  SQUARE 

.79666E*00 

■UMBER  OF  RESIDUALS 

93 

RESIDUAL  STAID ARD  ERROR 

89256E*00 

Figure  C17.  COWL83  Univariate  Parameter  Values 
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Figure  C19.  COWL83  Univariate  Histogram  of  Residuals 
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PfiEiWITENED  de*and  data 
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Figure  C20.  COWL83  Univariate  Power  Spectrum 
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SONUS!  OF  MODEL  1 

iiiitiiiiiiitititMiiiiiin  tut  iiiiiiiii  it  it  mini  tint  limit  iiiiiiiii  tiiimiiiiimiitmimtiimiiittiiiiiti 


DATA  -  Z  *  dcaand  data  ilv«8S  98  OBSEBVATIOIS 

DIFFEBEKIBG  01  Z  -  1  OF  OBOES  1 

•iittiiiiittiiiiiiiiiittiiiiiiititiiittiiitititiiititititiitimtiiitiiiiitititiintiiitmiiiiititimmiiiitmi 


PABAMETEB 

PABAMETEB 

PABAMETEB 

ESTIMATES 

98  PES  CEST 

IDMBES 

TYPE 

OBOES 

7  ALOE 

LOVES  LIMIT 

OPPEB  LIMIT 

IIIIIHIHIHHfflftllHIIIIIIIimiHtVtllHfllftlllltlllflflftHimHIIIIIIIIflllimfimifllllllltlflftffllfltltl 
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AOTOBEGSESSIVE  1 

1 

-.34146E-01 

*.23550E»00 

. 16721E*00 
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A0T0BE0BESSITE  1 
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. 1 1681E+00 

.54171E*00 
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IDTIMG  ATE  BASE  1 

2 
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-,12945E*00 

29591E*00 
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NDVIIG  A7EBAGE  1 
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-.49543E-0) 
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Figure  C21.  SLVE85  Univariate  Parameter  Values 
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Figure  C23.  SLVE85  Univariate  Histogram  of  Residuals 
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Figure  C24.  SLVE85  Univariate  Power  Spectrum 
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SUMMARY  OF  MODEL  1 
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DIFFERENCING  ON  2  -  1  OF  ORDER  I 

iilllHIlHmHIfHitKHmtitUHIItHKmHItltKWlimtiHitmHHtMWHtltfiHItmttWWmMttHII 

PARAMETER  PARAftTER  PARAMETER  ESTIMATED  95  PER'  CENT 

NUMBER  TYPE  ORDER  VALUE  .  LONER  LIMIT  UPPER  LIMIT 

1  MOVING  AVERAGE  1  0  382S8E+>'«C  -.57397E+00  -.19179E+00 

tHmtttHHHIItIHttHttHHHHmWttmttHWIHWHmttHHHmHUifHtmtWmHmtHIIWtHHW 

OTHER  INFORMATION  AND  RESULTS 

HittHWItlfHtmtHIHiHItHHiliiWHWWHIHWIHifHttitHilWWHtllltHtHHHtHWIHWtlHHIl 
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Figure  C25.  ASSY32  Univariate  Parameter  Values 
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Figure  C27.  ASSY32  Univariate  Histogram  of  Residuals 
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Figure  C28.  ASSY32  Univariate  Power  Spectrum 
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suttiary  of  model  1 
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DIFFERENCING  ON  Z  -  1  OF  ORDER  1 
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OTHER  INFORMATION  AND  .RESULTS 
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NUMBER  OF  RESIDUALS  95  RESIDUAL  STANDARD  ERROR'  .3I051E+01 


Figure  C29 .  ASSY33  Univariate  Parameter  Values 
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Figure  C31.  ASSY33  Univariate  Histogram  of  Kesiduals 
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SUMMARY  OF  MODEL 
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Figure  C33.  FORK35  Univariate  Parameter  Values 
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Figure  C35.  FORK35  Univariate  Histogram  of  Residuals 
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